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ABSTRACT 

Electromagnetic  interference  -from  power  lines  is  one  of 
■the  main  sources  of  man-made  interference  to  communications 
in  the  high  Frequency  (HF)  radio  band.  Two  types  of  radio 
interference  generated  by  power  lines  are  gap-type  noise 
caused  by  electric  discharges  across  line  hardware  and 
corona  noise  caused  by  the  partial  breakdown  of  the  air  due 
to  the  high  electric  fields  around  transmission  line 
conductors.  Using  original  data,  this  research  has  devel- 
oped a  parameter  based  model  of  gap-type  and  corona  noise 
that  allows  the  fundamental  noise  mechanisms  to  be 
mathematically  or  physically  simulated.  An  expression  for 
the  power  spectral  density  (PSD)  of  gap-type  noise  and 
corona  is  derived.  The  energy  detection  problem  is 
formulated,  and  using  analytical  results  based  on  the  Hall 
model  for  radio  noise,  a  robust  energy  detection  receiver  is 
developed.  Tests  of  this  receiver  using  actual  and 
simulated  data  are  described. 
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I .  INTRODUCTION 

A.  PURPOSE 

Observations  of  radio  interference  at  high  -frequency 
(HF)  receiver  sites  have  indicated  that  existing  models  of 
man-made  radio  noise  are  inadequate  to  describe  the  observed 
time-,  and  -Frequency-domain  behavior  oF  the  noise.  The  data 
also  shows  that  quite  often  a  specific  noise  source  pre- 
dominates in  a  particular  location.  These  observations  led 
to  the  objectives  of  this  dissertation: 

1.  to  develop  specific  models  for  certain  man-made 
noise  sources  in  the  HF  radio  band  and 

2.  to  apply  this  model  to  the  analysis  of  energy 
detection  receivers. 

B.  BACKGROUND 

Studies  conducted  over  the  past  few  years  have  indicated 
that  one  of  the  primary  sources  of  man-made  radio  noise  in 
the  HF  band  Are  alternating  current  (AC)  transmission  and 
distribution  lines  CRefs.  1,23.  Two  of  the  primary  sources 
of  power  line  noise  aure  gap  noise,  also  known  as  micro- 
sparking,  and  corona.  Gap  noise  is  caused  by  a  sparking 
process  between  hardware  points  on  utility  poles  and  corona 
is  caused  by  the  partial  breakdown  of  air  due  to  high 
electric  fields  around  high  voltage  conductors.  Both  of 
these  noise  types  are  non-Gaussian  noise  processes. 
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The  majority  of  the  models  of  man-made  and  atmospheric 
noise  to  date  have  assumed  the  high  amplitude  impulses 
driving  the  receiver  have  either  a  Poisson  arrival  rate  or  a 
variation  of  Poisson  arrivals.  This  assumption  has  allowed 
the  derivation  o-f  first  order  envelope  statistics:  the 
amplitude  probability  distribution  (APD)  -function  and  the 
density  function  of  the  phase  of  received  noise  for  various 
amplitude  distributions  of  the  driving  impulses.  Observa- 
tions of  man-made  noise  from  power  lines,  however,  has 
indicated  that  the  assumption  of  interpulse  independence 
(Poisson  arrival  times)  is  not  valid.  There  exists  a 
definite  time  domain  correlation  of  the  impulses  driving  the 
receiver.  This  is  due  to  the  underlying  deterministic 
mechanism  of  the  fundamental  frequency  of  the  power  line 
voltage.  Consequently,  this  research  has  concentrated  on 
the  statistics  of  the  impulse  arrivals  which  anr&  manifest  in 
the  autocorrelation  function  and  the  spectrum  of  the 
observed  noise  process.  Accordingly,  the  models  developed 
are  specific  rather  than  general  and  are  more  suited  to 
source  identification  by  spectrum  analysis  and  robustness 
evaluation  of  systems  rather  than  to  generic  specification 
of  optimum  receiver  structures. 

The  performance  of  energy  detection  receivers  for 
stochastic  signals  in  non-Gaussian  noise  has  received 
relatively  little  study.  Actual  implementation  of  algo- 
rithms for   signal  detection   has  been   based  on  heuristics. 


The  advent  of  digital  signal  processing  techniques  that 
allow  -For  complex  post-detection  algorithms  suggests  a 
careful  statistical  analysis  o-F  the  noise  performance  oF 
these  systems  may  help  develop  more  eF-Ficient  and  robust 
receivers . 

C.   CONTENT 

A  brief  outline  oF  the  rest  oF  the  thesis  will  now  be 
presented.  Chapter  II  develops  a  general  model  For  bandpass 
impulsive  phenomena  using  a  Filtered  impulse  model  and 
complex  envelope  theory.  Well  known  empirical  models  and 
physical  models  o-F  radio  noise  are  presented  and  their 
relationship  to  the  generalized  model  discussed.  Particular 
attention  is  paid  to  the  Hall  model  For  atmospheric  radio 
noise  CRe-F.  3D.  Some  useful  extensions  to  it  are  derived 
that  will  be  used  For  simulations  in  Chapter  V . 

In  Chapter  III,  Field  observations  From  sources  oF  gap 
noise  interference  are  presented  and  analyzed.  Three  cases 
oF  actual  interference  are  used  and  a  probabilistic  model 
based  on  a  statistical  analysis  of  the  data  and  the  filtered 
impulse  model  developed  in  Chapter  II  is  specified.  The 
power  spectral  density  (PSD)  of  gap  type  noise  as  predicted 
by  the  model  is  derived  and  compared  to  the  data. 

Using  an  actual  case  of  corona  noise  from  a  500  kv"  power 
line,  the  same  type  of  analysis  is  carried  out  for  corona 
noise.   The  corona  noise  model  is  based  on  the  same  filtered 
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impulse  framework  as  "the  gap  noise  model  of  Chapter  III 
however,  "the  specifications  differ  substantially. 

Using  the  Hall  model  -For  atmospheric  noise  discussed  in 
Chapter  II,  a  locally  optimum  receiver  -For  detection  of 
unknown  signals  in  HF  atmospheric  noise  is  derived.  A 
practical  modification  to  this  receiver  shows  it  to  be  an 
adaptive  limiter.  The  performance  of  this  receiver  is  then 
examined  in  simulated  man-made  noise  and  in  recorded  HF 
signal  and  noise  data. 

Appendix  A  presents  a  brief  description  of  the 
instrumentation  that  was  used  to  collect  much  of  the  data 
for  this  thesis.  The  dissertation  concludes  with  a  brief 
summary  of  results  and  some  suggestions  for  future  research. 


1  1 


II .  IMPULSIVE  NOISE  AND  RECEIVER  MODEL 

A.   INTRODUCTION 

In  "this  chapter  a  general  set  of  specifications  -For  an 
impulsive  noise  interference  process  will  be  developed  and 
the  results  of  previous  work  in  the  field  explained  in  terms 
of  the  specifications.  This  generalized  scenario  will  be 
used  in  later  chapters  to  specify  complicated  sources  of 
man-made  radio  interference.  In  attempting  to  describe 
atmospheric  radio  noise,  two  general  types  of  models  have 
been  developed:  empirical  models  designed  to  fit  first  order 
statistical  data,  and  physical  models  directly  related  to 
the  underlying  physical  mechanisms.  The  Hall  empirical 
model  for  atmospheric  radio  noise  has  been  shown  to  fit 
atmospheric  noise  data  very  well  and  will  be  used  in  this 
dissertation.  Some  extensions  to  the  Hall  model  will  be 
developed  in  this  chapter. 

8.   GENERALIZED  SCENARIO 

In  order  to  provide  a  framework  for  the  discussion  of 
the  impulsive  noise  models,  a  general  interference  scenario 
for  impulsive  noise  will  be  described  CRef.  4,53.  A  typical 
interference  scenario  consists  of  the  following  elements: 

1.  a  source  of  interference, 

2.  a  transmission  medium  to  the  receiver  and 
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3.  the  receiver  where  the  interference  manifests 
itself. 

For  this  study  the  noise  process  at  the  receiver  is 
modeled  as  the  sum  of  a  high  density  (in  time),  low  amp- 
litude Gaussian  component  and  a  low  density,  high  amplitude 
impulsive  term. 

These  elements  are  shown  below,  where  e'Ct)  in  Fig.  1  is 
the  impulsive  interference,  L'Coj)  is  the  frequency  response 
of  the  transmission  medium,  z  '  (t)  is  white  Gaussian  noise 
and  H'(u>)  is  the  combined  response  of  the  RF  and  IF  filters 
of  the  receiver. 


z  '  (t) 


e  '  (t) 


n  '  (t) 


Figure  1.   Interference  Scenario 
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The  source   of   impulsive   noise   is   specified   by   "the 

equation 

N(-t) 
e'(t)  =  E  e,6(t-t»)  (2-1) 

i  =  l 

where  et   is  the  real  amplitude  of  the  ith  pulse  and  N(t)  is 

a. unit   counting  process   that  generates   arrival  times,  the 

ti  's.   Using  this  representation,  the  impulsive  interference 

source  is   described  by   the   probability   density   -Function 

(PDF)  o-f   each  ei  ,   given  by   pt  (e)  ,  and  the  impulse  arrival 

times  generated  by  the  unit  counting  process,  N(t). 

In  order   to  simplify   later  analysis   the   interference 

scenario  will  be  expressed  in  terms  of  its  complex  envelope. 

The  interference   process  at   the  receiver  detector  is  n'(t) 

and  will   be  replaced   by  its   complex  envelope   equivalent, 

n(t) ,  where 

jo>©  t 
n'(t)  =  ReCn(t)e     ]  (2-2) 

and 

n(t)  =  nc (t)  +  jn, (t) .  (2-3) 

jU)0  t 

The   term    n(t)e  is   known   as   the   analytic  signal 

representation  of  n'(t),  and  n(t)  is  the  complex  envelope 
representation  of  n'(t).  The  reference  frequency  for  the 
complex  envelope  representation  is  cu0  .  When  n'(t)  is  a 
bandpass  process,  nc (t)  and  n« (t)  are  the  lowpass  inphase 
and  quadrature  modulation  components  of  n(t)  respect- 
ively.CRef.  6:p.  753 


14 


SpeciFying  a  -Filtered  impulse  model  in  complex  envelope 
■Form  requires  deriving  an  expression  For  the  driving  impulse 
Function  given  in  Eqn.  2-1.  Two  equivalent  approaches  may 
be  taken  and  they  are  both  outlined  below.  In  the  First 
approach  a  bandpass  impulse  may  be  postulated  where  the 
spectrum  oF  the  impulse  is  assumed  Flat  From  u)0  -  W  to  u)0  + 
W  where  u>0  is  the  reference  Frequency  and  W  is  an  arbitrary 
bandwidth.  For  this  case  the  analytic  signal  representation 
is  easily  -Found  in  the  Frequency  domain.  Using  the 
Frequency  domain  deFinition,  the  analytic  signal  associated 
with  e'(t)  in  Eqn.  2-1  is 


eA  <t)  =  (1/tt) 


CD 

jtut 
E'(u>)e    da).  (2-4) 


0 

Inserting  the   definition  o-F   the  bandpass   impulse   defined 

above  and  using  Eqn.  2-1 

N(t)       jtuo  (t-t,  ) 
eA (t)  =  E  (2/n)  e  sinCW (t-t»  ) ) / (t-t. ) .  (2-5) 

i  =  l 

The    complex    envelope    oF    e'(t)     is    defined    as 

-JCOo  t 

e(t)    =    eA (t)e  (2-6) 

and 

N(t)  -j8, 

e(t)     =    E     (2/n)    e  sin(W(t-t, ) ) / (t-tt  )  (2-7) 

i  =  l 

where  <D0ti   =  8»  .    As   W  becomes   large,  e(t)  will  asympto- 
tically approach  CRef.  71 


15 


N(t)   -j8, 
e(t)  =  E   2e      & (t-t, ) .  (2-8) 

i  =  l 

The  second  approach  CReF.  83  is  to  consider  the  impulse 
as  ideal  and  constant  in  the  -Frequency  domain.  The  analytic 
signal  representation  oF  e'(t)  is  then  directly  determined 
■From  the  time  domain  definition  oF  an  analytic  signal.  The 
analytic  signal  is  defined  as 

e* (t)  =  e'(t)  +  je'(t)  (2-9) 

where  e'(t)  is  the  Hilbert  transform  oF  e'(t).   The  analytic 

signal  representation   oF  Eqn.   2-1  using   the  above   deFin- 

itions  is 

N(t) 
eA (t)  =  E  6(t-tt)  +  j/(n(t-t»))  (2-10) 

i  =  l 

since  1/nt   is  the   Hilbert  transForm  oF  the  delta  Function. 

The  complex  envelope  oF  e'(t)  is  then 

N(t)         -j9i      -jtDot 
e(t)  =  E  6(t-t.)e     +  je      /(n(t-t4))       (2-11) 
i  =  l 

using  Eqn.  2-6. 

Although  apparently   different,  when  the  complex  impulse 

trains  described   by  Eqns.  2-8  and  2-11  are  convolved  with  a 

complex   Filter   impulse   response,   the   Filter   output   is 

identical  For   both  representations.   IF  h(t)  is  the  complex 

envelope  oF   the  Filter   impulse  response,   then  the   Filter 

output  with  the  input  given  by  Eqn.  2-8  is 
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N(t)  -j8, 

E   e     h(t-t,  ) 
1  =  1 

ihere  complex  convolution  is  defined  as 


02-12) 


h(t)     *    e(t)     =     ( 1/2) 


h(T  )e(t— DdT 


(2-13) 


The    filter    output    with    the     input    defined    by    Eqn.     2-11     is 


N(t)  -j9t  -juJot 

(1/2)    E    h(t-t, )e  +    je  hA  (t-t, ) 

i  =  l 


(2-14) 


where  the  convolution  definition  of  the  Hilbert  transform 
was  used  CRef.  6:p.  69].  Using  Eqn.  2-6,  the  same  result  as 
derived  in  Eqn.  2-12  is  then  found.  Due  to  its  simplicity 
and  for  ease  of  use  in  programming  simulations,  Eqn.  2-8 
will  be  used  to  represent  the  driving  impulse  function  for 
-the  rest  of  this  dissertation. 

One  further  point  to  mention  concerns  the  probability 
distribution  of  8t  in  Eqn.  2-8.  For  the  processes  we  con- 
sider, ti  will  have  random  arrival  times  and,  as  can  be  seen 
in  Eqn.  2-7,  9t  is  formed  by  multiplying  t,  times  the 
reference  frequency,  u)<>  .  Consequently,  9i  is  the  phase  of 
the  impulse  time  with  respect  to  the  reference  frequency. 
It  can  be  assumed  to  be  uniformly  distributed  over  0  to  2n 
when  the  tt  's  have  a  probabilistic  interarrival  distribution 
and  u)o  is  much  greater  than  the  inverse  of  the  interarrival 
times  CRef.  9:  pp .  279,  10]. 
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In  an  actual  interference  situation  the  received 
impulses  are  filtered  by  the  radiating  antenna,  spreading 
losses,  atmospheric  attenuation,  receiving  antenna  and 
cabling  losses.  These  terms  will  be  lumped  together  as  an 
equivalent  -Filter,  L  '  ( w )  .  Examples  of  this  -Filtering  term 
For  various  interference  sources  have  been  recorded  by  many 
researchers  CRef.  11,123.  The  Filter  characteristic  is 
specified  in  the  Frequency  domain  as  L'(oj)  since  that  is 
where  it  is  most  easily  observed. 

An  example  of  this  type  of  filtering  is  shown  in  Fig.  2 
where  the  frequency  spectrum  of  an  impulsive  noise  source 
due  to  a  gap  discharge  is  shown  from  0  to  200  MHz.  The 
straight  lines  in  the  3-Axis  view  are  due  to  stations  in  the 
high  frequency  (HF)  band  from  2  to  30  MHz  and  the  FM  band 
from  88  to  106  MHz.  The  continuous  envelope  seen  in  the 
upper  view  is  the  magnitude  of  the  frequency  response  of 
L '(a)).  The  solid  line  in  the  upper  picture  is  the  noise 
floor  of  the  spectrum  analyzer  at  -100  dBm . 

The  complex   envelope  of   the  impulse   response   of   the 

attenuation  term,  L '  ( U) )  ,  is  given  in  the  frequency  domain  by 

L(u)>  =  LP  CL  '  ( -(jj-ujo  )  +  L'"  (uj+Uo  ]  (2-15) 

where  LP   denotes   the   lowpass   part   of   the   quantity   in 

brackets . 

A  complex  Gaussian  component  z(t)  will  also  exist  in  the 
observed  noise   process  and   is  due   to  the   combination  of: 
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Figure  2.   Gap  Discharge  in  Frequency  Domain 
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1.  the  thermal  noise  in  the  receiver,  and 

2.  the  combined  sum  of  many  low  level  atmospheric  or 
man-made  impulsive  sources. 

This  term   will  be  modeled  by  independent  zero  mean  Gaussian 

noise  processes   -For  the  inphase  and  quadrature  terms  with  a 

power  spectral  density  (PSD)  in  each  component  o-P  N0/2.   The 

distinction  between   the  Gaussian   noise  term   and  impulsive 

noise  term   is  that   many  receiver  responses  overlap  for  the 

high  density,   low  amplitude   Gaussian  case,   satisfying  the 

condition  -for   the  Central  Limit  Theorem.   For  the  impulsive 

noise  the   receiver  response   to  each  impulse  is  discernible 

such  that   the  probability   o-P  more   than  one   or  two  pulses 

overlapping  is   negligible.    Common   sources   o-P   impulsive 

noise  are   man-made  noise  due  to  power  lines  and  atmospheric 

noise  due  to  lightning. 

The  time   domain  input  to  the  receiver  is  now  modeled  by 

the  expression 


N(t) 
n(t)  =  (h(t)/2)  *  Cz(t)  +  E  2et e 

i  =  l 


where 


jQi 


id)  ( t-ti  ) 


L(u>)du>] 


(2-16) 


h(t)  =  hc(t)  +  jhs(t). 
The  noise   power  bandwidth   o-P  H  ( u) )  will  be  designated  BEff. 
In  most   situations  L'(w)   will  effectively   be  constant   in 
comparison  to   the  narrowband   receiver  filter   H(u)).    This 
makes   L  ( cu )    a   constant    that    depends    only    on    oj0  . 
Additionally,  g(t)  will  be  defined  by 
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g(t)  =  (h(t/2))  *  z(t)  (2-17) 

and  is  "the  low  level  Gaussian  noise  a"b  the  detector.  The 
power  in  the  inphase  and  quadrature  components  of  the 
Gaussian  noise  will  be  N<,Bt»,/2.  The  complex  lowpass  noise 
process  can  now  be  written  as 


N(t)   -j8, 
n<t)  =  g(t)  +   E  a, e     h(t-tt) 
i  =  l 


(2-18) 


where 


a,  =  e,  L'  (cdo  )  •  (2-19) 

The  parameter  a,  is  then  a  scenario  dependent  parameter  that 
is  a  function  of  the  reference  frequency  and  the  impulse 
amplitude.  With  the  above  simplifications,  Fig.  1  can  now 
be  expressed  in  complex  envelope  form  as  shown  in  Fig.  3. 


a(t) 


z(t) 


>   n(t) 


Figure  3.   Complex  Envelope  of  Interference  Scenario 

The  terms  from  Fig.  1  have  been  replaced  by  their  equivalent 
complex  envelope  representations  and  the  prime  is  dropped  to 
indicate  this. 
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Some  further  s  1  mp  1  1 -Pi  cat  1  ons  will  now  be  considered 
based  on  the  characteristics  of  the  H  '  (  cj  )  ,  the  bandpass 
■Filter.  IF  the  impulse  response  h(t)  of  the  Filter  is  real, 
H'(cu)  has  conjugate  symmetry  about  the  origin  in  the  sense 
that  H'"(oj)  =  H '(-a)).  Furthermore,  lF  H'(o))  is  sym-metrical 
in  the  same  way  about  the  reference  Frequency,  u)<>  ,  then 
hs  (t)  will  equal  0. 

The  envelope  squared  oF  the  process  defined  in  Eqn.  2-18 

1  s 

2  2  2 

E  (t)  =  nc   (t)  +  n9   (t) .  (2-20) 

IF  the  Gaussian  noise  term  is  assumed  to  equal  0  then 

2      N(T)  N(T) 
E  (t)  =  E     E   a,  aj  hc  (t-t,  )  hc  (t-tj  )cos(9,  -0j  )  + 
1=1   j=l 

N(T)  N(T) 

E     E   at  aj  h9  (t-tt  )  h9  (t-tj  )cos(9,  -Gj  )  .  (2-21) 
1=1   j=l 

IF  N(T)   defines  a   low  density   counting  process   such  that 

there  is   negligible  overlap   between  subsequent   pulses  and 

h(t)  is   the  complex   envelope  of   the  impulse  response  of  a 

bandpass  symmetrical   filter,  then   the  envelope  of  n(t)  can 

be  simply  written  as 

N(t) 
E(t)  =  E  a,h(t-t, )  (2-22) 

i  =  l 

where  h(t)  is  equal  to  hc  (t) . 

At  this   point   of   the   analysis   two   avenues   may   be 

pursued.   If  there  is  a  time  dependency  between  pulses  (non- 
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exponential  ly  distributed  interarrival  times),  then  the 
correlation  function  and  spectrum  o-F  the  process  can  be 
examined.  If  the  amplitude  probability  distribution  (APD) 
of  the  envelope  is  the  desired  result,  then  the  char- 
acteristic function  of  the  envelope  is  most  useful.  The  APD 
function  is  commonly  used  in  radio  noise  research  and  is 
equal  to  one  minus  the  cumulative  distribution  function 
(CDF).  In  general  these  two  approaches  a^r&  mutually 
exclusive  since  the  assumption  of  time  dependency  between 
the  pulses  makes  the  envelope  APD  function  difficult  to 
calculate.  If  the  impulse  arrival  times  form  a  Poisson 
point  process  then  the  envelope  APD  may  be  calculated: 
however,  the  spectrum  is  constant.  The  remainder  of  this 
chapter  will  follow  the  envelope  APD  approach.  Results 
obtained  by  previous  researchers  will  be  developed  using  the 
general  noise  model.  Chapters  III  and  IV  will  explain  the 
other  avenue  and  look  at  the  spectrum  of  complicated  man- 
made  noise  processes. 

C.   EMPIRICAL  MODELS 

Empirical  models  of  atmospheric  noise  have  been 
developed  to  provide  a  mathematical  expression  for  the 
first-order  statistics  of  the  envelope  of  the  received 
waveform.  In  particular,  the  APD  has  been  emphasized.  This 
type  of  model  attempts  to  construct  a  mathematical  express- 
ion that   fits  observed   data  without  regard  for  the  physics 
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of  the  interference  scenario.  Empirical  models  have  the 
advantage  that  the  resulting  expressions  are  much  simpler 
than  those  obtained  -Prom  the  -Filtered  impulse  models.  One 
disadvantage  oF  this  class  o-P  model  is  that  since  only  the 
-First  order  statistics  oF  the  noise  are  considered  in 
developing  and  Fitting  the  model,  the  higher  order  sta- 
tistics may  not  match  the  data  well. 

Variations  oF  the  Rayleigh  distribution  have  been 
proposed  by  a  number  o-F  authors  to  -Fit  observed  atmospheric 
noise  data.  The  Rayleigh  PDF  is 

2 

-or 
p(r)  =  2are         r  2    0  (2-23) 

=  0  r  <  0 

and    the    APD    or    exceedance    probability 

2 
-ocr0 
Pr(r    >    r„  )     =    1     -    P(.r0)     =    e  r0     >    0  (2-24) 

=1  r,     <    0 

How  well  the  single  parameter  Rayleigh  distribution  Fits 

observed  atmospheric   noise  data  can  be  seen  in  Fig.  4.   The 

data  points  are  exceedance  probabilities  plotted  in  decibels 

(dB)  above  the  root-mean-squared  (RMS)  value  oF  the  received 

envelope.    The  data   was  measured  on  an  ARN-2  receiver  at  a 

center  Frequency   oF  10   MHz  at  Boulder,  Colorado  CRef.  133. 

The  high  probability,  low  amplitude  portions  oF  the  observed 

data  approaches  the  slope  oF  the  Rayleigh  distribution  curve 

For  a       equal  to   12.  (a  is  defined  in  Eqn.  2-24)   This  value 
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Figure  4.   Atmospheric  Noise  and  Rayleigh  Distribution 
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oP  a   was  chosen   to  Pit  the  high  probability  portion  oP  the 

data.    That  the   data  approaches  a  Ray  lei gh  distribution  is 

due  to   the  e-PPect   oP  the   many  small   overlapping  impulses 

occurring   at    low   amplitudes   whose   quadrature   Gaussian 

components  will   have  a   Rayleigh  envelope.   However,  at  the 

low  probability,   high   amplitude   part   oP   the   curve   the 

Rayleigh   distribution    predicts  -far  Pewer   values   than 

actually  occur.    This   is  the   distinguishing  Pactor  oP  the 

envelope   distribution    o-P   impulsive    noise:    that    its 

distribution  contains  a  higher  probability  o-P  high  amplitude 

terms  than   is  predicted   by  the  Rayleigh  distribution.   The 

second  curve   in  Fig.   4  where   a  is   equal   to   .5   is   the 

Rayleigh  distribution   with  the   same  power   as  the  observed 

data.    It  can   be  seen   that  this  curve  under  estimates  the 

high  amplitude   and  over   estimates  the  low  amplitude  values 

o-P  the   observed  distribution.   It  should  be  noted  that  this 

type  o-P   plot  emphasizes   the   low   probability   portion   o-P 

distribution  since   that  is   where  the   deviation   Prom   the 

expected  behavior  occurs. 

In  an   attempt  to   correct  the   poor  Pit  oP  the  Rayleigh 

distribution  at   high  amplitude   levels,  Likhter   [Ref.   14] 

proposed  a  combination  of  two  Rayleigh  distributions: 

2         2 
-ar0       -flr0 
Pr(r  >  r0)  =  (1  -  c)e      +  ce     .         (2-25) 

This  formula   has  been   shown  to  agree  poorly  with  practical 

results  CReP.  153. 
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Spaulding,  Roubique   and  Crichlow  CRe-P.  16]  combined  the 

Rayleigh  distribution   with  a   "power"  Rayleigh  distribution 

to  obtain   a  distribution  that  -Pit  very  well  for  atmospheric 

noise  over  a  wide  range  o-P  receiver  bandwidths.   This  APD  is 

2 
-ocr0 
Pr(r    >    r0  )     =    e  r0     <    fl 

2    1/s 
-(ar0     ) 
Pr(r    >    r0 >     =    e  r0     2    fl  (2-26) 

where  a,   fl  and   s  are   determined  -From  measured  statistical 

parameters  of  the  noise. 

Horner  and   Harwood  CRef .  17]  proposed  the  two  parameter 

log-normal  distribution  and  -Pound  it  gave  a  satisfactory  Pit 

to  radio   noise  data   in  the   VLF  band.    The  log-normal  was 

chosen  since   it  has  a  more  impulsive  tail  and  Pits  the  high 

amplitude,  low   probability  data   better  than   the   Rayleigh 

distribution.   The  PDF  for  the  log-normal  distribution  is 

2    2 
1  -dn(r)  -  a)  /2a 

p(r)  =  e  (2-27) 

1/2 
ra (2n ) 

Two  examples   o-P  the   APD  o-P  the  log-normal  distribution  ar& 

shown  in   Fig.  5  where  the  curve  with  a  equal  to  -.24  and  <r2 

equal  to   .24  is   the  log-normal   distribution  with  the  same 

mean  and  mean  square  parameters  as  the  curve  with  a  equal  to 

.5  in   Fig.  4.    Comparing   the  two   curves   the   log-normal 

distribution   provides    an   intrinsically    better   Pit   to 
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Figure    5.       Atmospheric    Noise    and    Log-Normal    Distribution 
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impulsive  noise  data  at  both  the  high  and  low  amplitude 
values.  When  the  log-normal  density  has  both  parameters 
chosen  to  match  the  mean  and  mean  square  values  o-f  the 
observed  data  an  excellent  -Fit  to  this  particular  data  is 
obtained.  This  is  seen  in  the  curve  with  a  equal  to  -1.4 
and  a*  equal  to   1.4. 

Another  model  that  has  been  successfully  used  to  -fit  the 
observed  APD  o-P  atmospheric  radio  noise  is  Hall's  general- 
ized "t"  distribution  CRe-f.  33.  This  model  is  unique  in 
that  it  is  not  o-f  the  -filtered  impulse  type;  however,  in 
contrast  to  the  strictly  empirical  models  it  does  postulate 
a  random  process. 

One  o-f  the  problems  o-f  Gaussian  models  o-f  atmospheric 
noise  is  that  the  dynamic  range  o-f  the  model  is  less  than 
the  observed  dynamic  range  o-f  the  noise.  To  achieve  a 
greater  dynamic  range  starting  with  a  Gaussian  process,  Hall 
proposed  a  model  which  considered  the  received  noise  to  be 
o-f  the  -form 

n(t)  =  a(t)s(t)  (2-28) 

where  a(t)  is  a  slowly  varying  stationary  random  process  and 
s(t)  is  an  independent,  narrowband  Gaussian  process.  An 
analytically  tractable  distribution  was  selected  -for  a(t) 
which  was  chosen  to  give  good  agreement  between  the  model 
and  measured  atmospheric  noise  data.  The  distribution  -for 
a  (t )  is: 
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2  2 
(m/2)  1      -(m/2a  a     ) 

p(a)  =  — — e 

m  m+1 

a    r(m/2)      |a|  (2-29) 

where  m   and  a       are    the  two  parameters  used  to  -Fit  the  model 

to  the  data.   The  distribution  oF  s(t)  is  Gaussian; 

2     2 

1  -(s  /2<r»   ) 
p(s)  =  e 

2  1/2 

(2n<r,  )  (2-30) 

Hall  calculated  the  density  -Function  oF  n(t)  to  be 

8-1 
T(9/2)      y  1 

p(n)  =  

1/2      2     2  8/2 
T((8-l)/2)    n        (n   +  y    )        (2-31) 

where 

1/2 
Y  =  m    (<^i  /<^5  » 

8  =  m  =  1  , 

and  r   is  the   gamma  Function.    For   the  case  where  at  =  a , 

Eqn.   2-31    is   the   density   Function   For   Student's   "t" 

distribution  which   is  the   basis  For  describing  the  density 

as  a  generalized  "t"  distribution. 

Hall  then   calculated  the  envelope  distribution  For  n(t) 

based  on  the  above  assumptions  and  obtained 

8-1  r 


p(r)  =  (8  -  1)  y 


2     2  ((8  +  D/2) 
(r   +  y  >  (2-32) 


with  the  phase  uniformly  distributed  between  0  and  2n 
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For  atmospheric  radio  noise,  values  of  9  in  the  range  2 
to  5  have  been  -Found  to  give  excellent  agreement  with  the 
data  CRef.  18,193.  (Models  with  integer  values  of  the 
parameter  0  will  be  abbreviated  to  Hal  12,  Hal  13,  etc.) 

Fig.  6  shows  the  Hal  13  -Fit  to  the  same  atmospheric  noise 
data  used  in  Figs.  3  and  4.  An  excellent  -Fit  to  this  data 
is  obtained.  The  Hall  model  has  some  disadvantages.  The 
higher  moments  only  exist  For  orders  greater  than  8-1  and 
the  parameters  9  and  y  must  be  determined  For  each 
interference  scenario  and  are  not  easily  related  to  the 
physical  source  of  the  interference. 

D.   FILTERED  IMPULSE  NOISE  MODELS 

The  Filtered  impulse  models  diFFer  Fundamentally  in 
concept  From  the  empirical  models  described  previously. 
This  class  oF  model  will  be  explained  in  terms  oF  the 
generalized  interference  scenario  presented  at  the  beginning 
of  the  chapter.  The  Filtered  impulse  models  have  the 
advantage  oF  being  based  on  the  underlying  physics  oF  the 
process  but  suffer  from  the  disadvantage  of  being 
analytically  complex.  Interestingly,  some  of  the  results 
obtained  from  the  filtered  impulse  models  have  retroactively 
justified  the  expressions  that  were  derived  from  the 
empirical  models. 

To  determine  the  APD  of  the  envelope  of  a  filtered 
impulse  process,   the  joint   characteristic  function   of  the 
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Figure  6.   Atmospheric  Noise  and  Hall  Distribution 


32 


inphase  and   quadrature  components  is  determined.   Since  the 

noise  is   narrowband  it   will  be   circularly  symmetrical  and 

the  2-D   Fourier  transform   will  become   a  Hankel   trans-Form 

CRef.  20,21].    This   characteristic  -Function  is  trans-Formed 

to  polar   coordinates,  then   inverted  to   determine  the  APD . 

Starting  with  the  complex  noise  process  described  in  Eqn.  2- 

18,  the   joint  characteristic   Function  of   the  inphase   and 

quadrature  components  is 

j  (u)i  nc  (t)  ■+•  oja  n3  (t)  ) 
0  Co),  ,oja  )  =  ECe  ]  .       (2-33) 

It  is   shown  in   Appendix  B  that  the  characteristic  Function 

oF  the  envelope  due  to  the  Gaussian  term  is 

2 
-N0B(cur  )/2 
09  (uv  )  =  e  (2-34) 

and  the   characteristic  Function   oF  the  envelope  due  to  the 

impulsive  term  is 


X 


p(a) 


[Jo  (uf  ahc  (t))  -  lDdTda 


(2-35) 


0,  (ojr  )  =  e 
where  p(a)  is  the  impulse  amplitude  density  Function,  X,  the 
Poisson  rate  Function,  hc  (t) ,  the  inphase  component  oF  the 
Filter  impulse  response  and  J0  is  the  ordinary  Bessel 
Function  oF  the  First  kind.  The  density  Function  oF  the 
envelope  is  then  the  inverse  Hankel  transForm  given  by 


p(r)  = 


ruv  Jo  (raj,  )0a  ( uv  )0i  ( u)r  )  du>r  . 


(2-36) 
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This  equation  For  -the  PDF  is  diFFicult  to  evaluate;  however, 
a  number  oF  researchers  have  examined  special  cases  to 
obtain  results. 

Furutsu  and  Ishida  CReF.  223  derived  an  equation  o-P  the 
same  -Form  as  Eqn.  2-36  with  no  Gaussian  term.  They  obtained 
many  approximate  results  by  considering  the  equation  -for 
speciFic  source  density  Functions.  IF  p(a)  is  exponential 
they  Found  that  the  resulting  envelope  density  was 
approximately  Rayleigh  For  small  amplitude  values.  IF  p(a) 
was  uniformly  distributed,  the  log-normal  density  was  Found 
to  be  a  good  approximation  over  most  oF  the  range  oF 
interest.  IF  a  strong  local  source  oF  interference  was 
present,  Furutsu  and  Ishida  showed  that  a  Function  involving 
the  confluent  hypergeometr ic  Function  was  a  good 
approximation.  They  also  considered  one  special  case  were 
the  unit  counting  Function  N(T)  was  modiFied  to  be  a  Poisson 
Poisson  branching  process  and  showed  that  For  low 
amplitudes  the  resulting  density  Function  was  approximated 
by  the  simple  Poisson  case. 

Giordano  CReF.  23D  evaluated  Eqn.  2-36  For  various 
distributions  oF  p(a)  which  he  determined  From  assumed 
propagation  laws  and  spatial  distributions  oF  impulse 
sources.  In  one  particular  case  he  adopted  the  Following 
assumptions : 

1 .  No  Gaussian  component 
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2.  Uni-Form  spatial  distribution  o-F  impulse  sources 

3.  Inverse  distance  law  -For  received  -Field  strength 

4.  Arbitrary  receiver  envelope  response 

and  showed  that  the  APD  o-P  the  envelope  -Function  was  approx- 
imated by 

K 

Pr(r  >  r0  )  = — — — — (2-37) 

2     2  1/2 
(r0   +  K  ) 

where 

T 

K  =  c  hCt)  dt  , 
0 
T  is  the  observation  interval,  h(t)  the  -Filter  response,  and 
c  is  a  scenario  dependent  constant.-  This  result  is  sig- 
nificant because  it  has  the  same  Form  as  the  APD  oF  the 
Hal  12  model  and  physically  justi-Fies  what  had  heretoFore 
been  an  empirical  model. 

Middleton  CReF.  24D  included  the  Gaussian  noise  term  in 
Eqn.  2-36  and  expanded  the  characteristic  -Function  without 
taking  the  expectation.  Then  by  inverting  the  char- 
acteristic function  term  by  term,  a  canonical  Form  oF  the 
envelope  density  Function  can  be  obtained.  The  exceedance 
probability  was  then  shown  to  be 

Pr(r    >    r<>  )    = 

m 
2  2<x>(-l)m  ma  ma 

-r0     CI     -    r0     E   A.     TCI     +  )»F,(1 ;     2;     r„  )  1 

m=0       m!  2  2 

e 

(2-38) 
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where  t  F,  is  a  confluent  hypergeometric  function,  and  o  and 
A.  are  -two  parameters  that  ar&  determined  by  the  source  dis- 
tributions and  the  propagation  law. 

E.   HALL  MODEL  EXTENSIONS 

For  later  work  in  this  dissertation  a  model  for 
atmospheric  noise  statistics  is  required.  This  section  will 
further  describe  the  Hall  model  as  applied  to  HF  atmospheric 
noise  and  will  derive  some  extensions  to  it  that  will  be 
used  to  digitally  simulate  atmospheric  radio  noise. 

The  Hall   model  for   the  envelope   statistics   of   radio 

noise  is 

9-1          r 
p(r)  =  (9  -  1)  y      


2     2  <(9-l)/2) 
(r   +  y  )  (2-39) 

and  was   chosen  for  a  number  of  reasons.   The  most  important 

aurez    the   Hall  model   is  analytically   tractable  and  shows  a 

good  fit  to  the  data.   It  can  be  obtained  for  certain  values 

of  the   parameter  8,   starting   from   the   physically   based 

filtered  impulse   models.    This  was   shown  by  both  Giordano 

CRef.  23D   and  Middleton  CRef.  24].   Additionally,  Schonhoff 

CRef.  18]   has  related   the  Hall   models  to   CCIR  Report  332 

CRef.  25],   the  standard  reference  of  first  order  statistics 

of  atmospheric  noise.   In  the  CCIR  report  the  APD  curves  are 

shown  as  a  function  of  Vo .   The  parameter,  V&  ,  is  a  commonly 

used  measure  of  the  impulsiveness  of  atmospheric  radio  noise 

and  is  defined  as 
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2  1/2 
ELr  ] 


VD  =  20  log  

ECr]  (2-40) 

where  r   is  the   amplitude  of   the  received   envelope.    For 

quadrature  Gaussian   noise  V0   =  1.05   dB  and   as  the   noise 

envelope  becomes  more  impulsive,  VD  will  increase. 

One  problem  in  relating  the  Hall  model  to  VD  is  that  the 

moments  of  the  Hall  model  only  exist  -For  orders  greater  than 

9-1.    This  means  that  for  the  Hall2  and  Hall3  models,  V0  is 

undefined  unless   the  model   is  modified.    Two  methods  have 

been  proposed   to  do   this.    Hall 's  method   CRef.   3]   adds 

another  term  to  the  envelope  density  function  and  multiplies 

by  a   negative  exponential   term  to   reduce  the   tail  of  the 

density  function   and  give  finite  moments.   The  advantage  of 

this  method   is  that   the  density   function  exists  from  0  to 

infinity;   however,    the   resulting    density   function   is 

complicated.    The  method   proposed  by   Schonhoff  CRef.   18] 

truncates  the   Hal  12  and   Hal  13  distributions   above   a   set 

level,  TP,   and  renormalizes  the  density  f unct ion  by  picking 

k  such  that 

TP 

k  p(r)dr  =  1 .  (2-41 ) 

0 

The  level,   TP  ,  is   chosen  to   set   a   desired   V0   for   the 

distribution.   The  constant  k  will  be  a  function  of  D  where 


2  1/2 
D  =  (1  +  (TP /y)  ) . 


(2-42) 
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The  existence  oP  the  -first  and  second  moments  -for  these  two 
distributions  is  then  assured.  Using  this  procedure 
Schonhoff  generated  a  -Family  oP  distributions  parameterized 
by  8  that  can  be  used  to  represent  a  wide  range  o-P  measured 
data.  By  matching  the  VD  ratio  o-P  the  data  to  the  dis- 
tribution, a  close  -Pit  to  the  APD  is  Pound. 

An  extension  to  the  basic  Hall  model  to  be  used  in  a 
later  chapter  is  the  density  Punction  oP'  a  constant 
amplitude  sinusoid  plus  the  assumed  noise.  This  density  is 
analogous  to  the  Rician  density  in  Gaussian  noise  theory. 
It  will  be  used  in  simulating  the  perPormance  oP  energy 
detection  receivers  Por.  the  signal  plus  noise  case.  It  was 
derived  in  CRe-P.  263  Por  the  Hal  12  model  and  it  is  derived 
Por  the  Hal  13  and  Hal  15  models  in  Appendix  C.  For  the  Hal  13 
model  the  density  Punction  is 

2         2         2  2 

tr    +    A    +    y    )     2y    rD 

p(r)     =    (2-43) 

4  22  22  22  242 

(r    -    2A    r       +    2A    y    2+    2y    r    2+    A    +    y    )  (D    -     1) 

where  A   is   the   signal   amplitude,  y       is   a   Hal  13   noise 

parameter  and   D  is   dePined  in   Eqn.  2-42.     Likewise   the 

density  Punction  Por  the  Hal  15  model  is 

4    2      2 
2y  r(b   +  2a  ) 

p(r)  =  (2-44) 

2     2  5/2 
(a   -  b  ) 
where 

2     2      2 
a    =    \r       +  A   +   y 
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and 

b  =  -2Ar. 
Fig.  7  is  an  example  oF  the  APD  For  a  sinusoid  plus  Hal  13 
noise  with  the  RMS  value  oF  the  noise  set  to  1  and  the 
signal  to  noise  ratio  at  -26,  -6  and  14  dB .  For  the  -26  dB 
SNR  case  the  APD  is  very  similar  to  the  Hal  13  noise  only 
density  shown  in  Fig.  6. 

Commonly  used  -Functions  relating  to  the  Hall  models  ar& 
shown  in  Table  I.  These  are  the  envelope  density  -Function, 
the  1st  and  2nd  moments  oF  the  envelope,  the  inverse  oF  the 
CDF,  the  quadrature  marginal  density  Function  and  the  Hall 
noise  plus  random  sinusoid  density.  The  parameter  D  in 
Table  I  is  deFined  by  Eqn.  2-42  and  the  Function  £C-3  is  the 
complete  elliptic  integral  oF  the  second  kind.  The  quad- 
rature components  can  be  derived  by  transForming  the 
envelope  density  to  rectangular  co-ordinates  and  integrating 
with  respect  to  one  component.  The  inverse  oF  the  envelope 
CDF  is  also  included  since  it  provides  an  easy  method  to 
generate  random  deviates. 

Also  useful   For  simulation   purposes  are   the  values  oF 
the  Hall   parameters,  8   and  y»   ^or  diFFerent   values  oF  v"D 
normalized  to  a  unit  root  mean  square  value.   The  truncation 
point  TP  was  Found  For  a  desired  VD  using  an  iterative  tech- 
nique.  These  parameters  are  shown  in  Table  II. 
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V,, 

Hall  No.   ECn  (t)  :   ECn(t)D 

TP 

Y 

D 

2.1 

5        3 

L         .785 

- 

1  .01 

- 

3 

4        3 

L         .708 

- 

.708 

- 

4 

3        3 

L         .631 

11  .66 

.420 

27.76 

5 

3        3 

L         .562 

25.90 

.364 

71  .05 

6 

3        1 

L         .501 

67.71 

.321 

211  .  1 

8 

2        3 

L         .398 

11  .77 

.085 

137.6 

10 

2        3 

L        .316 

17.08 

.058 

291  .0 

12 

2        1 

L         .251 

24.25 

.041 

587.  1 

14 

2        3 

L         .  199 

33.81 

.029 

1141 .9 

TABLE  II 

F .   SUMMARY 

Both  empirical  and  analytical  -Filtered  impulse  models 
have  been  successful  in  explaining  the  first  order  envelope 
statistics  of  bandpass  radio  noise  processes.  However,  in 
obtaining  these  results  the  assumption  has  to  be  made  that 
the  impulse  arrival  times  due  to  the  noise  process  form  a 
Poisson  point  process.  In  general  this  will  not  be  true  for 
man-made  impulsive  noise  and  -for  some  types  of  atmospheric 
noise  CRef.  27D .  In  most  cases  of  man-made  noise  the 
underlying  deterministic  mechanism  of  the  device  causing  the 
noise  will  modulate  or  regulate  the  impulse  distribution. 
The  Hall   model  for   atmospheric  noise   has  been  examined  in 
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detail  and  some  useful  extensions  to  it  have  been  derived 
that  will  allow  noise  and  signal  plus  noise  environments  to 
be  simulated. 
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Ill .  GAP  NOISE  MODEL 

A.  INTRODUCTION 

In  this  chapter  a  specific  source  of  man-made  radio 
noise  will  be  modeled  and  specified  in  terms  of  the  counting 
process,  N(t)  ,  and  the  probability  density  Function  (PDF)  oF 
the  impulses,  p(a)  .  These  specifications  were  introduced  in 
the  previous  chapter  to  define  the  impulsive  noise  model. 
The  noise  source  that  will  be  examined  is  gap  type 
discharges  and  is  commonly  Found  on  electric  power  dis- 
tribution and  transmission  lines.  Gap  noise  is  one  oF  the 
major  types  oF  interference  From  power  lines  and  is 
Frequently  observed  as  the  primary  interference  to  com- 
munication systems  operating  in  the  high  Frequency  (HF) 
radio  band. 

B.  GAP  NOISE  MECHANISM 

At  least  two  mechanisms  have  been  Found  by  which  a  gap 
discharge  process  can  occur  on  a  power  line.  The  resistance 
in  the  line  insulators  can  be  degraded  allowing  current  to 
Flow  through  the  insulator  base  and  creating  a  potential 
gradient  across  any  gaps  or  defects  in  the  insulator 
mounting  hardware.  A  second  way  in  which  a  potential  can  be 
created  across  an  air  gap  is  by  an  electro-static  coupling 
of  the  line  potential  to  isolated  hardware  on  the  pole.  In 
both  cases  the  potential  across  the  gap  is  discharged  by  the 
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voltage  breakdown  of  the  gap  and  the  resulting  rapid  current 
-Flow  or  spark.  This  process  generates  a  radio  Frequency 
(RF)  noise  impulse  with  spectral  components  extending  into 
the  hundreds  oF  MHz.  During  a  single  discharge,  the 
potential  across  the  gap  is  temporarily  diminished.  How- 
ever, while  the  Fundamental  60-Hz  waveform  is  still  greater 
than  a  threshold  voltage,  the  process  can  occur  again.  The 
spark  will  discharge  across  the  gap  .  repeatedly  until  the 
alternating  current  waveform  drops  below  the  breakdown 
threshold  potential.  CRef.  28:p.  78,293 

The  spark  discharge  and  recharging  oF  the  gap  potential 
indicates  that  this  type  oF  process  is  regenerative  and  can 
be  modeled  as  a  renewal  process  where  the  renewal  points  ar& 
associated  with  the  sparks.  One  Feature  oF  gap  noise  that 
complicates  the  modeling  is  the  120-Hz  on-oFF-on  modulation 
imposed  on  the  renewal  points  by  the  alternating  current 
(AC)  waveform.  One  way  to  account  For  this  effect  is  to 
consider  the  turn  on  times  as  another  renewal  process  driv- 
ing the  spark  discharge  process.  This  type  of  model  is 
known  as  a  branching  renewal  process  CRef.  30D  . 

The  noise  processes  that  were  used  to  develop  this  model 
were  short  term  stationary  in  that  the  statistics  of  the 
process  did  not  change  over  the  observation  interval .  Not 
all  sources  of  power  noise  are  stationary  even  over  a  short 
observation  interval.  Physical  effects  such  as  wind,  solar 
heating  and  varying  line  loads  can  act  to  make  certain  types 
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o-P  power  line  noise  highly  variable  -from  observation  to 
observation.  Although  these  types  of  noise  were  observed 
during  data  collection,  the  data  used  to  develop  the  model 
was  taken  -From  noise  sources  that  were  stationary  Por  the 
length  o-f  the  ten  minute  data  records. 

C.   OBSERVED  TIME-DOMAIN  CHARACTERISTICS 

The  gap  noise  sources  -for  this  study  were  observed  on 
utility  distribution  lines  in  the  vicinity  of  the  Naval 
Postgraduate  School  and  were  chosen  to  illustrate  parameters 
of  the  noise  model.  Fig.  8  is  a  typical  time-domain  observ- 
ation of  a  gap  discharge  process  observed  at  3  MHz  and 
envelope  demodulated  with  a  10-kHz  Gaussian  bandpass  Pi  Iter. 
The  important  characteristics  to  note  aret 

1.  the  process  has  an  on-oPP-on  modulation  at  a  120- 
Hz  rate  related  to  60-Hz  waveform  of  the  power 

1 ine ,  and 

2.  the  pulse  groups  that  result  Prom  the  modulation 
have  a  probabilistic  number  oP  impulses  occurring 
in  each  group  and  a  probabilistic  interarrival 
time  between  pulses  in  a  group. 

In  the   Pirst  pulse   group,  9   pulses  occur  with  varying 

amplitudes  and   interpulse  arrival   times.    In   the   second 

pulse  group  starting  approximately  8.33  msec  later,  8  pulses 

occur,  again  with  varying  amplitudes  and  interarrival  times. 

One  group   oP  pulses   is  associated  with  either  the  positive 

or  negative   polarity  oP   the  line   voltage  wavePorm  and  the 

other  group   with  the   opposite  polarity.   Ident i Picat ion  oP 
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Figure    8.       Envelope    o-F    Demodulated    Gap    Noise 
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the  polari-ty   is  impossible   without  physically  locating  the 
source . 

The  average  amplitude  of  the  impulses  will  be  one  of  the 
parameters  -for  the  model.  An  estimate  of  this  parameter  can 
be  made  -From  the  above  presentation.  It  is  important  to 
note  that  -For  this  highly  impulsive  type  o-F  noise  that  the 
observed  amplitude  is  a  Function  oF  the  shape  and  width  o-F 
the  effective  bandpass  Filter  and  the  detector  char- 
acteristics. 

At  this  point  it  should  be  noted  that  the  model  can  be 
approached  at  two  levels  oF  complexity.  In  the  simple 
version  the  data  -From  both  the  positive  and  negative  phases 
is  considered  as  a  whole  and  averaged  over  the  two  phases  to 
determine  a  set  o-F  average  parameters.  In  the  more  complex, 
but  more  accurate  model,  two- sets  oF  model  parameters  are 
determined  -  one  For  each  phase.  For  the  remainder  oF  this 
work  the  more  complex  case  will  be  considered.  The  results 
For  the  simple  case  are  presented  in  Moose  and  O'Dwyer  [Ref. 
313  and  can  easily  be  determined  From  the  complex  case  by 
setting  the  parameters  equal  For  both  phases  and  simplifying 
the  resulting  equations. 

In  order  to  better  characterize  the  interpulse  arrival 
times,  which  will  be  used  to  determine  two  model  parameters, 
the  rising  raster  capability  described  in  the  instru- 
mentation section  was  used  to  generate  the  display  of  a 
different  gap   noise  source   which  is   shown  in   Fig.  9.   In 
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Figure    9.       Multiple    Scans    o-F    Gap    Noise    I rrter-Perence 
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■this  picture   the  amplitude   oF  each   individual  record   was 

reduced  and  thresholded  so  that  all  that  essentially  remains 

is  the   time  oF   impulse  arrival   data  -For   28  time  records. 

The  time   base  o-f  the  display  was  intentionally  synchronized 

to  the   power   line   to   -Facilitate   taking   data   and   this 

accounts   -For   the   regularity   oF   the   pulse   groups   From 

observation  to   observation.    In  this  view  the  interarrival 

times  For   56  pulse  groups  (28  groups  o-F  one  polarity  and  28 

o-F  the   opposite  polarity)   can  be  determined  along  with  the 

number  o-F   pulses  in  each  oF  the  56  successive  pulse  groups. 

The  average  number  o-F  pulses  per  group  considered  separately 

For  each   phase  will   also  be  used  as  estimates  o-F  two  model 

parameters. 

Fig.  10   is  a   histogram   o-F   the   distribution   oF   the 

interarrival  times   between  the   observed  impulses   For   the 

same  noise   process  shown   in  Fig.   9.    To  use  this  data  to 

generate  parameters   -For  the  noise  model,  the  histogram  will 

be  approximated  by  a  continuous  density  Function.   The  gamma 

density  Function  CReF.  32: p.  18] 

r-1  -(Xt) 
p(t)  =  X(Xt)    e       /(r-1)!  (3-1) 

was  chosen   because  it   showed  a   close  Fit  to  the  data,  and 

its  characteristic   Function,  which   will  be   used  in   later 

derivations,  was   particularly  simple.    The   characteristic 

Function  oF  the  gamma  density  is 

-r 
0(ju>)  =  (1  -  jo)X)    .  (3-2) 
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The  r   and  X   parameters  are  simply  interpreted  in  terms 

of  an   underlying  Poisson   process  as   the  time   to  the   rth 

point  of   a  Poisson   process  of   intensity  X.    If   jj  is  the 

random  variable   assigned  to  the  observed  interarrival  times 

then  the  unbiased  estimate  o-P  the  mean  o-P  p  is 

_  N 

est(p)  =  (1/N)  E  p,  (3-3) 

i  =  l 

and  the  unbiased  estimate  o-P  the  variance  is 

2       N    2         _  2 
estCa,,  )  =  -C  E  pt  -  N  est(p)  >/(N-l).  (3-4) 

i  =  l 

The  unknown  parameters  oP  the  desired  gamma  density  Punction 

anr&    estimated   using  the   method  oP  moments  CRe-P.  33 :  p .  250] 

using  the  Pol  lowing  relations; 

_  2        2  2 

r  =  est(p)  /est(aM  )    and   X  =  est (p ) /est (a„  ).     (3-5) 

Using  this   approach,  the   gamma  density   Punction   used   to 

approximate  the   interarrival  time  histogram  is  also  plotted 

on  Fig.   10.    Figs.  11  and  12  are    based  on  the  same  data  as 

Fig.  10.   However,  they  ar&    separated  according  to  different 

phases.   It  can  be  seen  Prom  the  r  and  X  parameters  in  Table 

III  that  there  is  a  signiPicant  difference  between  phases  of 

the  same   gap  noise   process  shown   in   the   histogram   even 

though  the   mean  value  of  each  data  set  was  almost  the  same. 

The  data   in  Fig.   11  was   more  clustered   than  the   Fig.  12 

data.    This  results  in  a  higher  r  and  lower  X  parameter  for 
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the  same  mean  value.  Figs.  13  and  14  are  additional  sets  of 
gap  noise  data  that  were  observed  during  the  research. 

Table  III   summarizes  the   r  and   X  parameters  -Pitted  to 
the  data  described  above. 

r  X 

Case  1  33.2  37,000 

Case  1  Phase  A  63.4  71,400 

Case  1  Phase  B  23.3  25,900 

Case  2  44.3  74,600 

Case  3  14.1  19,200 

TABLE  III 

A  -Final  model  parameter  will  be  called  T*  and  is  a  constant 
delay  or  offset  o-P  every  second  pulse  group  that  exists  with 
respect  to  the  fundamental  voltage  waveform. 

In  summary   the  inputs  to  the  model  based  on  time-domain 
data  ar&i 

1.  An  estimate  o-P  the  average  amplitude  o-P  all 
observed  pulses. 

2.  An  estimate  of  the  average  number  o-P  pulses  per 
group  Por  each  phase. 

3.  An  estimate  o-P  the  mean  oP  the  interarrival  time 
o-P  the  pulses  Por  each  phase. 

4.  An  estimate  o-P  the  variance  o-P  the  interarrival 
time  between  pulses  -Por  each  phase. 

5.  An  estimate  o-P  the  delay  in  the  start  time  oP  the 
pulse  groups  o-P  one  phase  with  respect  to  the 
fundamental  voltage  waveform. 
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The  average  amplitude  of  pulses  -Prom  di  Pferent  gap  noise 
sources  is  highly  variable,  ranging  -Prom  the  instrumentation 
noise  floor  to  higher  than  any  observable  signal  in  the  HF 
band.  The  average  number  of  pulses  per  group  has  been 
observed  to  vary  Prom  1  to  greater  than  20.  The  mean  of  the 
interarrival  times  ranges  Prom  .1  to  1  ms  and  the 
coefficient  o-P  variation  (standard  deviation  divided  by  the 
mean)  oP  the  interarrival  times  has  ranged  Prom  approx- 
imately .1  to  .3  Por  the  gap  noise  processes  we  have 
observed . 

D.   MODEL  DEVELOPMENT 

The  principal  goal  oP  this  chapter  is  to  describe  a 
noise  model  that,  with  suitable  choice  oP  parameters,  pro- 
vides an  adequate  simulation  oP  the  actual  physical  noise 
mechanism  o-P  gap  noise.  In  this  section  the  time-domain 
data  and  parameters  derived  in  the  previous  section  will  be 
integrated  into  the  Pi ltered  impulse  model  developed  in 
Chapter  II.  Recalling  Prom  Chapter  II,  a  Pi ltered  impulse 
noise  process  is  given  in  complex  envelope  notation  by 

N(t)   -j8, 
n(t)  =  g(t)  +  E  a,  e     h(t-t»)  (3-6) 

i  =  l 

where  p(a)   is  the   amplitude  distribution   oP  the  impulses, 

N(t)  is   a  unit   counting  process  whose  statistics  determine 

the  impulse   arrival  times,  8  is  a  random  variable  uniPormly 
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d  1  str  1  buted  between  0  and  2n  and  h(t)  is  -the  impulse 
response  of  the  receiver  filter. 

In  order  to  physically  justify  a  filtered  impulse  model, 
the  impulse  duration  must  be  small  compared  to  the  inverse 
bandwidth  of  the  receiver  filter.  This  condition  is  easily 
met  in  the  case  of  gap  discharges.  Laboratory  analysis  o-P 
temporal  characteristics  of  gap  discharges  for  various 
geometries  shows  that  the  impulse  durations  range  from  10  to 
100's  of  nano-seconds  CRef.  34],  Therefore,  for  filter 
bandwidths  up  to  1  MHz  the  output  noise  process  will  only  be 
a  function  of  the  incident  time  of  the  impulse  and  the 
filter  response,  not  the  waveform  o-P    the  impulse. 

The  specification  of  the  statistics  of  N(t) ,  the  count- 
ing process  that  drives  the  model  is  then  crucial  to 
obtaining  an  accurate  representation  of  the  physical  noise 
process.  Based  on  the  electrical  characteristics  of  a  gap 
discharge  and  the  observed  time-domain  characteristics,  the 
following  assumptions  are    made; 

1.  A  primary  series  of  event  times  separated  by  an 
interval  T0 /2   exists.   T0  is  the  fundamental 
period  of  the  power  line  voltage  waveform.   The 
distribution  of  the  time  to  the  first  primary 
event  is  uniform  over  0  to  T0  . 

2.  A  subsidiary  process  commences  at  the  primary 
event  times  l .e.  T0 /2  ,  To   ,  3T0 /2   . . .   The 
subsidiary  process  is  a  renewal  process  that 
continues  for  Nl  or  N2  renewals  alternating 
between  Nl  and  N2  renewals  at  successive  primary 
poi nts . 
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3.  The  interarrival  times  to  the  -first  and  subsequent 
points  of  the  subsidiary  renewal  process  for  each 
phase  are  all  independent  and  identically 
distributed  (IID)  within  a  pulse.   Each  phase  has 
its  own  set  o-f  parameters. 

4.  The  weighting  or  amplitude  distribution  of  all 
impulses  in  the  subsidiary  process  is  a  constant. 

5.  A  constant  offset  or  delay  designated  by  TD  may 
exist  at  every  second  primary  renewal  point. 

The  -first  assumption  is  supported  by  the  periodicity  in 
the  data  that  is  related  to  the  -fundamental  frequency  of  the 
power  line.  The  large  majority  o-f  observed  gap  noise  sources 
had  pulse  groups  on  both  the  positive  and  negative  phases  of 
the  fundamental  waveform.  This  accounts  for  the  T0 /2  per- 
iodicity . 

The  second  assumption  addresses  the  differences  between 
the  sparking  phenomenon  on  the  positive  and  negative  phases 
o-f  the  fundamental  .  It  is  frequently  observed  that  the 
ECNpot  3  is  different  than  ECNNE#3  where  N* o  ■  and  NN E a  are 
the  random  number  of  discharges  in  the  pulse  groups.  This 
can  be  due  to  two  effects;  an  asymmetrical  gap  geometry  and 
the  fundamental  physical  difference  between  the  sparking 
mechanisms  for  positive  to  ground  and  negative  to  ground 
sparks.  Two  deterministic  integer  constants;  Nl  and  N2 
nearest  to  E[Np0« ]  and  ECN*K«]  respectively  are  assumed  as 
the  approximation  to  the  variable  number  of  sparks  per  half 
cycle  of  the  fundamental  waveform. 
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The  -third  assumption  is  based  on  the  -Fact  that  the 
"inception  of  gap  discharges  in  natural  air  and  the 
development  of  electron  avalanche  are  -fundamentally 
probabilistic  processes  that  depend  on  atmospheric  pressure, 
humidity,  presence  of  natural  ions,  electrode  surface  and  so 
on"  [Ref.  34].  In  view  of  the  above  statement  and 
considering  the  empirical  data,  the  justification  for  devel- 
oping the  gap  discharge  as  a  probabilistic  process  is  well 
founded.  The  assumption  of  independence  from  discharge  to 
discharge  is  not  as  well  justified.  Effects  such  as 
electrode  heating  after  the  initial  discharge  in  a  cycle 
could  act  to  make  the  average  interarrival  time  vary  from 
pulse  to  pulse  within  a  pulse  group.  However,  to  develop  a 
tractable  model  the  assumption  of  independence  between 
arrival  times  is  made.  The  gamma  density  function  deter- 
mined by  the  estimate  of  the  mean  and  variance  of  the 
interarrival  times  for  each  phase  is  used  to  define  the  PDF 
of  the  interarrival  times. 

The  fourth  assumption  concerns  the  amplitude  of  the 
impulses.  In  most  cases  the  amplitude  of  the  impulses 
within  a  pulse  was  nearly  constant.  In  some  cases  there  was 
a  difference  in  the  average  amplitude  from  negative  phase  to 
positive  phase,  however  to  keep  the  model  tractable  the 
amplitude  of  all  impulses  is  assumed  constant. 

The  fifth  assumption  incorporates  an  effect  that  will 
explain  one   of  the   features  observed  in  the  power  spectral 
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density  (PSD)  of  the  observed  process.  The  physical 
rationale  behind  this  assumption  is  that  there  may  exist  a 
difference  in  the  threshold  between  the  positive  and 
negative  phases,  thus  consistently  delaying  or  off sett  i  ng 
the  start  of  the  renewal  process  on  one  of  the  phases. 

The  above  assumptions  place  this  model  in  a  class  o-P 
processes  known  as  branching  processes  with  the  main  process 
being  a  degenerate  renewal  process  with  an  interarrival  PDF 
of 

p(t)  =  6 (t-T0 /2)  (3-7) 

The  subsidiary   process  is   a  renewal   process  o-P   Nl  or   N2 
points.    Consequently,  the  impulse  source  Por  the  gap  noise 
model  is  completely  specified  by  the  following  parameters: 
a  -  Amplitude  of  the  impulses, 

N1,N2  -  Number  of  impulses  in  the  subsidiary  process, 
ut  ,Ma  _  Mean  value  of  impulse  interarrival  time, 
Var(pt  ) ,Var(u8 )  -  Variance  of  interarrival  times, 
To     -  Offset  on  one  phase. 

To  complete  an  interference  scenario  N0Beff/2  and  h(t)  must 

also  be  specified.   The  expression  for  the  branching  renewal 

process  is  then 

M(T)   N.       -j8.„ 
n(t)  =  g(t)  +  E     E  a. „  e      h(t-t.„-T.  )        (3-8) 
m=l   n=l 

Fig.  15a, b   shows  sample   realizations  of   the  envelope, 

E(t),  of  the  noise  process  defined  above. 


62 


E(t) 


E(t> 


M(-t) 


ILL 


k  K  K 


o  -^ 


■»  t    (a 


-»  t 


TD 

^ 

1     1 

— — — ^ 

(b 


»  t    (c 


Figure  15.   Envelope  and  Counting  Function  oF  n(t) 

The  Gaussian  term  was  set  to  0  and  h(t)  was  approximated  as 
a  decaying  exponential.  Note  that  in  comparison  to  an 
actual  noise  process  illustrated  in  Fig.  8  the  amplitude  oF 
the  impulses  is  constant  with  a  -Fixed  number  oF  impulses  on 
the  positive  and  negative  polarity  oF  the  line  voltage 
wave-Form.  Fig.  15c  is  the  primary  counting  process,  M(T), 
driving  the  subsidiary  process  shown  in  Fig.  15b.  Fig.  15b 
includes  a  non-zero  TD  and  its  relationship  to  M(T)  is  shown 
in  Fig.  15c . 

Fig.  16  shows  a  realization  oF  the  envelope  oF  this 
process  For  Nl  and  N2  =  3,  T  =  30  ms  and  gamma  density 
Function  parameters  rl  and  r2  =  32  and  XI  and  \2  =  37,000. 
Random  interarrival   times  with   a  gamma   distribution   were 
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approximated  by  summing  exponentially  distributed  random 
variables  CRe-P.  353.  The  impulse  amplitudes  were  normalized 
to  one  and  a  background  level  oP  quadrature  Gaussian  noise 
at  No B/2  =  .0001  was  added  which  simulates  either  receiver 
noise  or  high  density,  low  amplitude  impulsive  noise. 

E.   COMPARISON  OF  CALCULATED  AND  OBSERVED  DATA 

One  partial  description  oP  a  noise  process  is  the 
amplitude  probability  distribution  (APD)  o-F  the  envelope. 
This  description  has  been  used  extensively  in  the  analysis 
o-P  noise  and  in  deriving  optimum  receivers.  For  the  N(t) 
specified  in  the  previous  section  this  calculation  would  be 
diPPicult.  A  second  noise  process  descriptor  is  the 
spectrum  o-P  the  noise  process.  This  descriptor  is  suited  to 
our  model  where  the  structure  oP  the  noise  is  contained  in 
the  counting  process  that  drives  the  impulse  generation.  In 
addition  to  receiver  noise  perPormance  evaluation,  spectral 
analysis  can  also  be  used  Por  noise  source  ident 1 Picat 1  on 
and  isolation  CReP.  363. 

In  order  to  determine  the  spectrum  oP  the  envelope 
squared  oP  the  branching  renewal  process  dePined  by  Eqn. 
3-8,  the  spectrum  o-P  the  envelope  squared  of  a  non-branching 
renewal  process  is  determined  in  Appendix  D.  The  results 
obtained  in  Appendix  D  ar&  then  used  to  derive  the  spectrum 
o-P  the  the  branching  renewal  process  which  is  postulated  to 
model  gap   noise.      Starting   with   the   branching   renewal 
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process  given   in  Eqn.   3-8,  and   assuming  g(t)   is  equal  to 

zero,  it   is  shown   in  Appendix   E  that   the  average   oP  the 

Bartlett  estimate  oP  the  PSD  is: 

2    5 
S(u),T)  =  |H«(u))|  /T  Z    Term  <n)  (3-9) 

n=l 

where  n   is  an   index  to   the  5   Terms  de-Pined  below.   These 

five  Terms  involved  in  the  summation  in  Eqn.  3-9  are  de-Pined 

as  Pol  lows : 

Term  1  -  The  inter-pulse  group  summations  between  pulses  in 
pulse  groups  with  Nl  pulses 

Term  1  =  M2  5<u>,Nl)  (3-9a) 

Term  2  -  The  inter-pulse  group  summations  between  pulses  in 
pulse  groups  with  N2  pulses 

Term  2  =  M2  S(<d,N2)  (3-9b) 

Term  3   -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  Nl  pulses 

cos(M2    2nu>/u)0  )     -    1 

Term    3    =    -C  -    M2    >    ft<a>,Nl)  (3-9c) 

cos  (2na)/u)0  )     -    1 

Term  4  -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  N2  pulses 

cos(M2  2nu)/u)0  )  -  1 

Term  4  =  -C -  M2  >  n<w,N2)      (3-9d) 

cos  (2nu)/(i)0  )  -  1 

Term  5   -  The  intra-pulse  group  summations  between  pulses  in 

groups  with   Nl  pulses  per  group  to  pulses  in  groups  with  N2 

pulses  and  vice  versa 
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cos(M2    2no)/(D0  )     -     1 
Term    5    =    2cos  (  (u)(T0  /2-T0  )  )  -C — — — — —  }\y(0t   ,  Nl  ,0a  ,  N2) 

cos  (2n(i)/u0  )    "-    1 

(3-9e) 

where  0,   and  0a  are  the  characteristic  -Functions  associated 

with  the   interarrival  times  -For  each  phase.   The  parameters 

Nl  and   N2  are       the  number   oF  renewals  associated  with  each 

phase  and  5,  Q ,  and  \y  are  de-Pined  below: 

2  N 

N  -  N|0Cj(d)  |  0(j(D)  (1  -  0Cju))  ) 

1(0, N)  = -   2  Re-C __ yf 

2  2 

|l   -  0(jO))|  (1   -  0(jU)))  (3-10) 


N+l  2 

|0(  jU))   -  0  (  j(l))      | 

0(0, N)  =  , 


2 
1  -  0(  jco)  |  (3-11) 


and 


Nl+1  N2+1 

(0t  (j(D)  -  0i  (j<u)     )  (0a  (-jtD)  -  0a  (-jw)     ) 

H/(0»  ,N1  ,0a  ,N2)  = . 

(1-0,  ( j(D)  )   (1   -  0a  (-JU))  ) 

(3-12) 
As  a  test  of  the  model  the  PSD  o-F  the  gap  noise  process 
shown  in  Fig.  9  was  determined  using  a  spectrum  analyzer  and 
also  analytically  determined  using  only  the  model  parameters 
From  the  time-domain  data  and  the  methodology  oF  the 
previous  section.  Fig.  17  is  a  computer  plot  oF  the 
analytic  estimate  o-F  the  PSD  which  is  compared  with  the 
observed  PSD  computed  on  a  Wavetek  UA500A  spectrum  analyzer 
shown  in  Fig.  18. 
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Figure  18.   Observed  Estimate  o-P  Power  Spectral  Density  (PSD) 
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Figure  17.   Predicted  Power  Spectral  Density  (PSD) 


Comparing  the  two  PSD's,  it  is  seen  that  the  analytical 
expression  correctly  predicted  all  of  the  significant 
■Features  seen  in  the  observed  PSD. 

The  components  of  the  spectrum  due  to  the  different 
terms  is  of  interest  and  offers  some  insight  into  the 
origins  o-P  the  features  of  the  observed  spectrum.  The  dom- 
inant feature  of  the  observed  spectrum  are  the  periodic  120 
Hz  spectral  lines.  At  800  Hz  the  spectral  lines  change,  and 
appear  as  odd  harmonics  of  60  Hz,  however,  still  remaining 
at  120  Hz  intervals.  Terms  3-5  Ars  responsible  for  the 
these  spectral  line  features.  Each  of  these  terms  is  a  comb 
function,  determined  by  the  fundamental  frequency  and 
observation  interval,   multiplied  by   an  envelope  determined 
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by  the  irvterarr i val  "time  distribution.  Figs.  19  and  20  are 
Terms  3  and  4  -  respectively  and  clearly  show  the  phase  to 
phase  difference.  The  more  peaked  behavior  of  Fig.  19  at 
1100  Hz  is  due  directly  to  the  more  clustered  interarrival 
time  behavior  on  one  phase.  This  effect  was  observed  in  the 
time  domain  in  Fig.  11.  When  the  terms  shown  in  Figs.  19 
and  20  are  added  together,  they  form  60  Hz  harmonic  spectral 
lines  and  when  Term  5,  shown  in  Fig.  21,  is  added,  positive 
reinforcement  occurs  at  120  Hz  and  negative  at  60  Hz 
harmonics.  This  produces  the  observed  120  Hz  harmonic 
spectral  lines.  When  the  cancellation  is  not  complete  due 
to  differences  in  the  characteristic  function  from  phase  to 
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Figure  21 . 


PSD  Term  due  to  Phase  A  to  Phase  B  Intra-Pulse  Group 
Cor re  1 at l ons 
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phase,  small   60  Hz   components  can  be  present;.   This  effect 
is  seen  in  both  the  observed  and  predicted  spectra. 

The  observed  spectra  changeover  -Prom  120  Hz  harmonics  to 
odd  60  Hz  harmonics  at  800  Hz  is  also  correctly  predicted  by 
the  model.  This  effect  is  due  to  the  factor,  TD  ,  appearing 
in  Term  5  which  causes  the  odd  60  Hz  harmonics  to  be 
positively  rein-forced  and  the  120  Hz  harmonics  to  be 
canceled  out  above  800  Hz.  The  -Pinal  effect  shown  in  the 
model  and  seen  in  the  observed  spectra  is  the  presence  o-P 
continuous  spectral  components.  The  continuous  spectra  due 
to  Terms  1  and  2  is  shown  in  Figs.  22  and  23  and  also 
differs  from  phase  to  phases  The  continuous  spectra  peaks 
markedly  at  1100  Hz  and  this  effect  is  apparent  in  the 
observed  spectra  as  a  rise  in  the  noise  floor.  The  con- 
tinuous spectral  term  arises  when  Term  1  is  added  to  Term  3 
and  Term  2  to  Term  4.  The  peak  value  at  1100  Hz  in  Terms  1 
and  2  is  greater  than  the  negative  values  in  Terms  3  and  4 
at  the  same  frequency  and  when  added  together  causes  a 
continuous  rise  in  the  spectral  floor.  The  reason  for  this 
is  that  the  l  to  i  +  1  and  l-l  interarrival  times  within  a 
pulse  group  are  more  correlated  than  the  corresponding 
interarrival  times  between  different  pulse  groups.  The  more 
correlated  arrivals  appear  in  Terms  1  and  2  and  cause  the 
spectral  peak  to  be  larger. 
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F.   SUMMARY 

A  nine  parameter  model  -For  a  single  source  of  gap  noise 
interference  was  developed  and  shown  to  accurately  predict 
the  PSD  of  narrow  bandwidth  envelope  demodulated  gap  noise 
interference.  If  less  accuracy  is  desired,  phase  to  phase 
differences  can  be  neglected  and  the  resulting  model  has 
only  five  parameters.  The  model  was  specified  in  terms  of 
an  impulse  driving  function,  and  techniques  to  estimate  the 
model  parameters  from  time  domain  data  were  described. 
Since  the  model  specifies  a  complex  driving  impulse  function 
it  can  be  used  as  a  simulated  interference  source  for 
arbitrary  filters  and  receivers. 


74 


IV.  CORONA  NOISE  MODEL 

A.  INTRODUCTION 

In  addition  to  gap  noise  which  was  discussed  in  Chapter 
III,  another  major  source  oF  power  line  radio  noise  is 
corona  noise.  Corona  is  diFFerent  in  many  respects  -From  gap 
noise,  and  the  model  used  to  describe  it  will  take  a 
diFFerent  approach.  In  contrast  to  gap  noise,  which  is 
observed  on  both  electric  power  transmission  and  dis- 
tribution lines,  corona  noise  is  observed  only  on  power 
transmission  lines  and  generally  has  a  Fundamental  Frequency 
oF  180  Hz  corresponding  to  the  three  phases  oF  a  60  Hz 
alternating  current  (AC)  system.  Corona  noise  is  caused  by 
a  partial  breakdown  oF  air  surrounding  a  conductor  which  is 
at  a  high  potential.  Consequently,  the  impulsive  structure 
oF  corona  noise  is  not  as  well  defined  as  the  highly 
impulsive  structure  oF  gap  noise.  Corona  noise  From  AC 
sources  appears  somewhat  like  amplitude  modulated  Gaussian 
noise. 

B.  CORONA  NOISE  THEORY 

When  ar>  increasing  potential  difference  is  applied 
between  two  electrodes,  a  breakdown  voltage  is  reached  that 
is  characterized  by  the  transition  oF  air  From  a  poor 
electrical  conductor  to  a  good  conductor.  IF  the  Field  in 
the  gap   between  the   two  electrodes   exceeds   the   electric 
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strength  of  air,  which  is  about  30  kV/cm  -For  gaps  over  1  cm 
in  length,  prior  to  the  spark  breakdown  voltage,  then  corona 
will  exist.  Sharp  points  and  small  radius  wires  are  avoided 
whenever  possible  on  power  lines,  and  corona  generally  does 
not  become  a  significant  source  of  radio  interference  until 
voltages  exceed  about  70  kv"  CRef.  28:p.  843. 

As  the  voltage  on  an  AC  line  is  increased,  corona  will 
appear  and  be  associated  with  either  the  positive,  negative 
or  possibly  both  half-cycles  of  the  fundamental  waveform. 
For  an  aged,  clean  transmission  line  it  appears  that  the 
threshold  voltage  for  negative  corona  is  less  than  that  of 
positive  corona.  If  the  transmission  line  has  any  unusual 
characteristics,  the  threshold  voltages  of  positive  and 
negative  corona  become  difficult  to  categorize  CRef.  37]. 
In  particular,  it  has  been  observed  that  positive  streamer 
corona,  which  is  more  disruptive  of  communications  fre- 
quencies than  negative  corona,  predominates  on  newly 
installed  lines  and  during  periods  of  precipitation.  For 
modeling  purposes,  the  important  characteristics  of  corona 
noise  A\r&  not   a   precise   description   of   the   physical 

mechanisms  but  the  relationship  between  the  observed  inter- 
ference waveform  and  the  line  voltage. 

The  procedure  most  commonly  used  for  the  prediction  of 
interference  due  to  corona  noise  gives  the  noise  power  as  a 
function  of  the  characteristics  of  the  interference 
scenario.    The  random   properties  of   the  corona   noise  are 
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averaged  out  CRef.  383.  These  characteristics  would 
commonly  include  distance  -Prom  the  line,  the  line  voltage 
and  a  weather  correction  term.  Each  of  these  terms  is 
experimentally  verified  and  related  to  some  physical  process 
in  the  generation  of  corona  noise.  While  this  procedure 
gives  an  estimate  o-P  the  noise  power  in  a  specified 
bandwidth,  it  does  not  incorporate  the  probabilistic  nature 
or  time-domain  periodicity  observed  in  the  noise. 

The  approach  taken  in  this  chapter  will  continue  with 
the  filtered  impulse  theme  of  the  dissertation.  For  corona 
noise,  the  counting  process,  N(T),  will  be  modeled  as  a 
periodically  modulated  stream  of  impulses.  This  approach 
will  have  the  advantage  of  being  accurately  able  to  simulate 
the  time-  and  frequency-domain  behavior  of  corona  noise. 
However,  as  with  the  gap  noise  model  it  will  suffer  from  the 
disadvantage  that  for  each  scenario,  a  set  of  parameters 
must  be  experimentally  determined. 

C.   TIME  DOMAIN  CHARACTERISTICS 

The  time-domain  data  for  corona  noise  was  obtained  using 
the  instrumented  van  described  in  Appendix  A.  The  antenna 
was  a  3-meter  whip  mounted  on  the  van  (which  was  parked 
underneath  a  3-phase  500-kV  transmission  line).  The  weather 
conditions  at  the  time  of  measurement  were  clear  and  sunny. 
Records  of  the  corona  noise  were  obtained  using  the  envelope 
detected  output   of  the   HP-141T  described  in  Appendix  A  and 
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recording  the  data  on  a  20-kHz  bandwidth,  85-dB  dynamic 
range  digital  audio  system.  The  recording  process  allowed 
subsequent  analysis  in  the  laboratory. 

A  37-ms  portion  of  a  data  record  is  shown  in  Fig.  24. 
The  data  was  taken  at  a  300-kHz  center  frequency  and  10-kHz 
intermediate  -Frequency  (IF).  This  particular  picture  was 
obtained  by  playing  back  the  tape,  re-digitizing  the  data 
and  displaying  it  with  a  signal -analysi s  software  package. 
One  of  the  important  Features  of  the  time-domain  record  in 
Fig.  24  is  the  periodicity.  The  -Fundamental  Frequency  oF 
180  Hz  can  be  determined  From  a  periodic  rise  in  the  mean 
value  o-F  the  noise  and  From  a  periodic  increase  in  high- 
amplitude  impulses.  The  Fundamental  noise  Frequency  oF  180 
Hz  on  a  three  phase  line  indicates  that  the  dominate  noise 
is  generated  on  one  halF-cycle  o-F  the  Fundamental  60-Hz 
waveform  on  each  oF  the  three  phases  oF  the  line.  This  is 
consistent  with  low-level  negative  corona  predicted  For  aged 
lines  in  clear  weather. 

The  periodogram  spectrum  oF  the  data  record  described 
above  is  shown  in  Fig.  25.  The  time-domain  observations  oF 
the  180-Hz  periodicity  ar&  confirmed  and  a  360-Hz  harmonic 
of  the  fundamental  is  also  present.  At  higher  frequencies, 
no  harmonics  ar&  observed  and  the  spectrum  becomes  white. 
This  is  in  distinct  contrast  to  the  gap  noise  described  in 
Chapter  III  where  the  spectrum  of  gap  noise  has  harmonics 
well  into   the  kHz   range.   The  periodicity  can  also  be  seen 

78 


1    _      l._ 


l_ J.. 


r  Q 

O 

i     CVJ 


00 


£* 

o? 

°-UJ 

in 


Q     cd 


I- 

CD 


T 1 

cd    a 


a 


3CniINDUW    C3ZI"Ih)WdON 


u 

a 

h- 

5 


a 
a  uj 
a  n 

CE 

z: 
a 
o 


u 

Ld 

in 

Q 

I 

UJ 
CD 

LD 
CD 


li 

a 
o 

£K 
Ld 

a. 


CL 

cc 


— i     © 


01 

£ 
0 

a* 

£ 


Hi 
0 

z 

c 

0 

L 
0 

u 


<3- 

CVJ 

5. 
3 
01 


79 


L I. ... 


.1    .  ...  L 


,  !...„ 


X 1 


r   CS 
CD 


Si 


h- 
CE 

13  a 

LYUJ 
Lie 

a.  en 

CO 
hhZ 

2£ 

O 

£c 
oG 

a:  2: 

LU 

a: 


03 

— — •      p— 
Z 

— M     h-« 

o 

"Z3!    0. 


N 

X 


a 
+■ 

UJ 
CO 

• 

GO 

i 
li 


JO 
T5 


! 


a    a    g    a 
— «    oj    en 

i       1       : 


£3 


t       -  r~ 

LO       CD 

i       1 


(3 
i 


~f|       a 

,,  a  N 


CD 


cr 


2       2 
i        2 


on 


Z 

o 

H 
3 
_1 
O 

cn 
ui 
a 

> 

CJ 
2Z 

LU 
ID 
O 
U 
QC 
U. 


03 


O 

Q. 

ii 
N 

X 

Q 
+ 
UJ 
C3 

a 
a 


II     0J 


u 

"Z 

UJ 

a 

UJ 
K 
Lu 

QC 

LU 

O 

a 
u 
o 


N 

I 


0 

c 
<u 

3 

u 

L 
LL 


— i      o 


£ 
0 

>» 

0 

c 
11 

3 

a 

s. 

LL 


ll 
Ul 

■-* 
0 

z 

c 

0 

5. 
0 

u 


ID 

ai 

s. 

3 
0> 


qp  a3ziidwaoN 


80 


in  Fig.  26  (which  is  the  same  data  recorded  on  the  3-Axis 
display).  The  contribution  oP  each  phase  oP  the  power  line 
to  the  total  corona  noise  process  can  be  clearly  seen  in  the 
lower  picture  in  Fig.  26. 

Another  important  -Feature  o-P  the  data  Por  modeling 
purposes  is  that  the  receiver  response  to  an  individual 
impulse  overlaps  with  the  response  oP  preceding  impulses. 
Although  the  noise  is  still  impulsive,  the  simpliPying 
assumption  oP  non-overlapping  pulses  is  not  valid.  It  is 
seen  in  Fig.  24  that  the  amplitude  oP  the  impulses  is  random 
and  there  is  no  Pine-grain  interarrival  time  structure  other 
than  the  periodicity  discussed  above.  Both  oP  these 
observations  are  in  direct  contrast  to  gap  noise  where  the 
amplitude  Prom  impulse  to  impulse  is  almost  constant,  and 
there  is  a  dePinite  impulse  interarrival  time  structure. 
These  observations  provide  a  method  to  distinguish  between 
gap  noise  and  corona  noise. 

An  interesting  aspect  oP  noise  Prom  transmission  lines 
was  noted  while  taking  measurements  Por  the  corona  model. 
It  was  diP-Picult  to  obtain  a  recording  oP  power-line  noise 
in  the  vicinity  oP  high  voltage  transmission  lines  due  only 
to  corona  noise.  In  many  cases,  a  gap  noise  component  was 
superimposed  on  the  corona  noise.  This  can  be  seen  by 
comparing  Figs.  26  and  27.  Fig.  27  is  a  view  oP  both  gap 
and  corona  noise.  The  corona  noise  amplitude  peaks  at  the  2 
uV  level,   and  the   amplitude  oP   the  gap-noise   impulses  is 

81 


—   3 


Ampl it 


0 


time  -  m< 


50 


Figure  26.   Corona  Noise  as  shown  on  3-Axis  Display 
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Figure  27.   Gap  and  Corona  Noise 
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about;  3.5  pV.  The  two  lower  pictures  are  oP  identical  data 
with  the  threshold  varied  to  emphasize  each  noise  type.  The 
lower  picture  has  the  threshold  set  at  about  3  pV 
eliminating  the  180-Hz  corona  noise  and  showing  the  gap 
noise  process  occurring  at  a  120-Hz  rate.  This  corresponds 
to  a  gap  noise  source  sparking  on  the  positive  and  negative 
halF-cycle  oF  a  single  phase  oF  a  3  phase  line.  When  the 
threshold  is  lowered  to  about  1  pv",  the  180-Hz  corona  noise 
dominates  as  seen  in  the  middle  picture  oF  Fig.  27.  This 
phenomenon  is  due  to  corona  discharge  on  the  negative  hal-P- 
cycle  oP  each  phase  oP  a  3-phase  system.  Without  the  abil- 
ity to  examine  the  time-  and  Prequency-domai  n  behavior  o-P 
the  noise  using  a  scanning  analyzer  and  3-Axis  display,  it 
would  be  very  diPPicult  to  di-PPerent iate  between  gap  and 
corona  noise  and  make  accurate  measurements  oP  either. 

Fig.  28  illustrates  corona  noise  with  diPPerent  phase- 
to-phase  characteristics.  In  this  example  the  noise  power 
due  to  the  corona  is  unequal  Prom  phase  to  phase.  The  noise 
occurring  on  the  dominant  phase  is  approximately  twice  the 
peak  voltage  oP  the  phase  with  the  lowest  amplitude. 

D.   CORONA  NOISE  MODEL 

Based  on  the  above  observations  o-P  the  demodulated 
envelope  oP  corona  noise,  a  speciPic  Filtered  impulse  model 
based  on  Eqn.  2-18  will  be  developed.  This  model  is  similar 
to   models   used   in   synthetic   hydrology   to   account   For 
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Figure  28.   Corona  Noise  -  Uneven  Phase  to  Phase 
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seasonal  variations   in  stream-Flow  data  CRef.  393.   Assuming 

no  Gaussian   noise  component,   a   lowpass   complex   -Filtered 

noise  process  can  be  written 

N(t)   -jGt 
n(t)  =   E  at e     6 (t-tt  ) .  (4-1) 

i  =  l 

Since  corona  noise  is  a  high  density  process,  the  envelope 
cannot  be  simply  defined  as  in  Eqn.  2-22  since  the  response 
of  the  -Pi  Iter  to  adjacent  pulses  will  overlap.  To  cir- 
cumvent this  dif-Piculty  it  will  be  assumed  that  the  envelope 
o-P  n(t)  can  be  modeled  by 

N(t) 
e(t)  =  E  atg(t-tt)  (4-2) 

i  =  l 

where  both  e(t)  and  g(t)  are  real,  positive  and  g(t)  depends 

on  h(t)  and  the  process  intensity. 

The  assumptions  to  dePine  the  model  are  as  Pol  lows; 

1.  The  rate  parameter  Por  the  counting  function  N(t) 
is  periodically  modulated  at  frequency  oj0  and  will 
be  approximated  by  a  finite  number  of  terms  in  a 
Fourier  series. 

2.  The  amplitude  distribution  of  successive  pulses 
given  by  the  a,  's  are    independent  and  identically 
distributed  with  density  function  p(a). 

3.  The  start  time  of  the  periodic  modulation  is 
uniformly  distributed  over  the  period  T0  where  T0 

=   1  /<!)<>  . 

Assumption  one  is  motivated  by  the  periodicity  observed  in 
the  data.  In  most  cases  the  fundamental  period,  u)0  will  be 
180  Hz.  If  the  noise  is  more  dominant  on  one  phase  of  the 
line,  as   shown  in   Fig.  28,   then  60   Hz  components   of  the 
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noise   will    exist   and   for   a   more   accurate   model   the 

-fundamental  period   would  be   60  Hz.    The  rate  parameter  is 

given  by 

NT 
X(t)  =  E  a„cos(o)0nt)  (4-3) 

n=0 

where  NT   is  the   order   at   which   the   Fourier   series   is 

truncated . 

The  second   assumption  is   arbitrary   in   terms   of   the 

distribution  -Function   chosen   -For   the   impulses.      It   is 

obvious  that   random  impulse   amplitudes  are  being  generated 

by  the   cumulative   effect   oF   the   corona   noise   sources. 

However,  it   would  be   difficult  to  accurately  fit  a  density 

function  to   the  data   due  to  the  time  varying  nature  of  the 

process.    The  independence   of  the  amplitude  bursts  is  also 

conjecture  but   facilitates  the   development  of   a  tractable 

model.   The  Rayleigh  density  function  is 

2 
-a  /2fl 
p(a)  =  (a/fl)e  a  >  0 

=0  a  <  0  (4-4) 

and  was  picked  to  model  the  random  amplitudes.  It  is 
analytically  tractable  and  can  be  justified  by  assuming  the 
quadrature  components  of  the  impulse  amplitudes  are  in- 
dependent Gaussian  events. 

The  final  assumption  is  made  because  in  most 
interference  scenarios,  the  power  line  noise  process  will 
not  be   synchronized  to   the  process   being  interfered  with. 
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Thus  the  phase  will  be  effectively  randomized  over  one 
fundamental  period.  In  some  of  the  observations  of  corona 
noise  (see  Figs.  26-28)  this  assumption  was  clearly  violated 
and  the  instrumentation  was  intentionally  synchronized  to 
the  power  line  for  clarity  of  presentation. 

The  model  could  be  further  generalized  by  allowing  the 
amplitude  density  function  p(a)  to  be  a  periodic  function  of 
time,  p(a,t)  independent  of  the  rate  parameter.  This  added 
complexity  would  make  the  model  more  accurate;  however,  more 
parameters  would  have  to  be  estimated. 

The  above  assumptions  classify  this  process  as  a 
compound  non-homogeneous  filtered  Poisson  process.  A  con- 
venient tool  to  study  these  processes  is  the  cumulant 
generating  function  CRef.  41 :p.  117] 

se(t) 
uy(s)  =  ln<ECe       D>  (4-5) 

where  In  is  the  natural  logarithm.  The  joint  cumulant  gen- 
erating function  is 

Si  e  (ti  )  +  sa  e  ( ta  ) 
H/(s,,sa)  =  ln<  ECe  ]>.      (4-6) 

For  the  noise  process  defined  by  the  first  two  assumptions, 
it  is  shown  in  Appendix  F  that 


;ah  (t-a  ) 


iy(s)  = 


p(a) 


X(cc)  C  e 


-  lDdada 


(4-7) 


and 
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l|/(Si    ,Sa   )      = 


p(a) 


Stah(tt-ot)    +    saah(t8-o) 
\(a)     C    e  -     lDdada 


-a>  -<d  (4-8) 

From    the       proper-ties    O'f    the    cumulant    generating    -Functions    it 
can    be    shown    that 

ECe(t)  1    =    \p'  (0)  ,  (4-9) 

VarCe(t)]    =    ^''(0)  (4-10) 

and 


xi    iy(0  ,0) 


Cov (ti , t8 )    = 


■OSt  TDSa 


(4-1 1 ) 


It    is    straight -Forward    to    show    that 


ECe(t)D    =    ECaU 


00 

NT 

E    a„  cos  (o)o  na  )  h  (t-a  )  da  , 
n=l 


(4-12) 


VarCe(t) ]    =    ECa    1 


"       NT  2 

E    a„  cos  (o)o  na  )  h     (t-a)da 
J     n=l 


(4-13) 


and 


CovCti  ,ta  3    =    ECa    1 


r      NT 

E    a„  cos(u)0  na)  h(ti  -a  )  h  (t8 -a)da  . 
n=l 


(4-14) 
The  autocorrelation  Function  can  be  Found  in  terms  o-F  the 
above  Functions  as 

R(t, ,ta )  =  Cov(t, ,ta )  +  ECe(t, ) 3ECe(ta ) ] .     (4-15) 
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The  above  equations,  prior  to  the  application  of 
assumption  3,  define  a  wide  sense  cyclostat ionary  process 
CRef.  40D  where 

ECe(t  +  mTo)]  =  ECe(t)]  (4-16) 

and 

R(t,  +  mTo  ,ta  +  mTo  )  =  R(t,  ,t«  )  .  (4-17) 

The  statistics  of  the  process  are  invariant  to  a  shift  of 
the  time  origin  by  integral  multiples  of  the  period,  T0  . 
This  characterization  is  particularly  appropriate  when  the 
process  is  observed  in  synchronization  with  the  -fundamental 
power  line  -Frequency. 

When  assumption  3  is  applied,  the  process  start  times 
are  randomly  shifted  over  one  period,  the  process  becomes 
wide  sense  stationary  and  the  mean  and  autocorrelation 
become 

To 


ECn(t)D  = 


ECn(t)]  dt 


(4-18) 


0 


and 


R(t)  = 


To 
R(t+T ,t)dt. 


(4-19) 


0 


The  power  spectral  density  (PSD)  of  the  process  is 


S((j)  = 


_       -jU)T 

R(T)e     di 


(4-20) 
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The  -Following  algorithm  is  an  example  of  how  the  noise 
process  specified  above  can  be  digitally  generated.  For  the 
algorithm  specified  below;  N0  is  the  number  of  samples  in 
one  period  T0  ,  <nt  >  is  the  real  sequence  to  be  generated,  i 
is  the  index  variable  for  the  sequence  and  j  and  k  are  dummy 
index  variables. 
1  .  i  =  0 

2.  n(i )  =  0 

3.  Generate  a  Rayleigh  random  variable,  R,  with 
parameter  fl  and  set  n(i)  =  n(x)  +  R 

NT 

4.  Set  intensity  X  =  E  a„  cos(u)0  ni  /N0  ) 

n=0 

5.  Generate  a  Poisson  random  variable,  P,  and  set  j  = 
trunc(P/\  +  .5) 

6.  Set  n(k)  =  0  for  k  »  i  +  1  to  j 

7.  i  =  j 

8.  Go  to  Step  3 

The  above  algorithm  will  continuously  generate  a  modulated 
Poisson  impulse  process  that  may  be  digitally  -Filtered 
concurrent  with  the  impulse  process  generation  or  subsequent 
to  the  generation  of  a  complete  time  record. 

Two  additional  points  should  be  noted.  A  number  oP 
impulse  amplitudes  could  stack  up  on  each  other  iP  the 
truncation  operation  returns  a  zero  value.  This  corresponds 
to  the  physical  case  where  multiple  impulses  occur  and  are 
unresolved  within  one  sampling  period.  BePore  a  record  o-P 
noise  samples  generated  by  the  above  process  is  used  a 
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truncated  uniform  random  variable  between  0  and  N0,  the  sampled 
period,  should   be  subtracted  -From  i.   This  ensures  that  the 
start  time  is  uniformly  distributed  over  one  period. 

E.   PARAMETER  ESTIMATION 

Having  specified  the  model,  the  task  of  determining  the 
parameters  -For  a  given  interference  scenario  remains.  The 
parameters  to  be  determined  are  the  <an>,  the  Fourier  Series 
coefficients  of  the  periodically  varying  rate  function  and 
the  A  parameter  in  the  probability  density  function  of 
impulse  amplitudes.  As  an  example,  the  parameters  for  the 
noise  process  described  in  the  time  domain  observations 
Section  will  be  determined. 

One  method  of  estimating  these  parameters  is  by 
considering  the  predicted  mean  and  variance  of  the  process 
as  given  by  Eqns.  4-12  and  4-13  and  matching  these 
parameters  to  the  data.  In  order  to  do  this  analysis,  the 
corona  source  has  to  be  observed  as  a  cyclostationary 
process.  This  requires  synchronizing  the  data  collection 
instrumentation  to  the  line  frequency  so  that  the  same  point 
in  each  cycle  can  be  examined.  This  is  difficult  to  do  for 
field  measurements. 

Another  option  for  determining  the  model  parameters  is 
the  predicted  stationary  mean,  variance  and  PSD  or  auto- 
correlation function.  The  PSD  as  defined  in  Eqn.  4-20  is 
stationary  so   no  line   synchronization  is   required  and 


92 


instrumentation  is  readily  available  to  make  a  periodogram 

estimate  oF  the  PSD. 

For  the  measurements  taken  in  this  research,  the  lowpass 

impulse  response  oF  the  -Filter  was  approximately 

2 
-at 
h(t)  =  (a/n)e  t  2    0 

=0  t  <  0.        (4-21) 

With  the  above  deFinition  the  -Filter  has  unity  gain  at  u) 
equal  to  0.  The  effective  Filter  response  g(t)  was  approx- 
imated by  h(t) . 

In  Appendix  G  it  is  shown  that 

2  2 
-o)o  n 


a„  cos(u>o  nt)       (4-22) 


1/2  NT        4a 
ECe(t)D  *  (fln/2)     E   e 

n=0 

1/2 
where  (nfl/2)      is  the  First  moment  oF  p(a),  the  probability 

density   Function  oF  the  impulse  amplitudes.   The  stationary 

mean  is 


1/2 
ECe(t) 1    •  (fln/2)     a0 . 

The  covariance  oF  e(t)  is  shown  to  be 


(4-23) 


-a(t, -t, )  /2 


Cov(t,  ,tt  )  =   (a/n)2fle 


a>  2 

r   NT  -2au 

E  a„  cos(o)o  n(u+(tt -»-tt  ) /2)e     du  . 
n=0 
t, -t. |  (4-24) 
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Using   Eqns.   4-15   and   4-19   and   substituting   the   above 
results,  the  stationary  autocorrelation  -Function  is 

2 

-OCT  /2  1/2  2 

R(t)  =  2fi  Ce        a«  (a/8n)    +  (l/4)a0 


2  2 
-u)o  n 


NT        4a    2 
+  (1/8)   E   e         a0  cosCuo  nT) ] . 
n=l 

(4-25) 

Taking  the  Fourier  transform,  the  resulting  PSD  is 

2 
1/2  -w    /2a  2 

S(uj)     =    2fl     C(2)  a0     e  +     (n     /2)a«     6(u>) 


NT  4a  2 

+     (n     /8)       £       e  a„     6  ( cu    ±     ncj0  )  . 

n=l 

(4-26) 

The  terms   in  the   PSD  in   order  are;   (1)  a  white  noise  term 

multiplied  by  the  -Frequency  response  oF  the  Gaussian  -Filter, 

(2)  an   impulse  at   cj  =  0  -From  the  mean  value  o-F  the  process 

and  (3)   periodic  components   at  harmonics   o-F  180  Hz  due  to 

the  modulation  o-F  the  driving  -Function. 

To  estimate   the  model   parameters,  the   sample  mean  and 

mean  square  estimates  were  used  in  conjunction  with  Eqns.  4- 

23  and  4-24  to  solve  For  a0  and  13.   The  periodogram  estimate 

o-F  the   spectrum  was   used  to  determine  the  magnitude  oF  the 

at     and   a2  co-e-F-Fi  c  l  ents .    The  Fourier  series  was  truncated 
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at  NT  =  2  since  no  higher  order  terms  were  observed  in  Fig. 
25.  For  the  data  taken  in  Fig.  24  and  25  the  parameter  a  in 
Eqn.  4-21  was  7  .  1  1  x  1  08  .  Table  IV  lists  the  parameters  of  the 
assumed  model. 


Coefficient  o-P  Gaussian  Filter 


7. 1 lxl08 


Amplitude  Density  Parameter 


fl 


1 .58xl02 ' 


Fourier  Co-efficients  of  Rate 
Funct l on 


a0  22,300 

a,  5,700 

aa  2,500 


TABLE     IV 

These  parameters  were  inserted  into  the  algorithm  described 
previously  for  generating  corona  noise,  and  five  cycles  of 
the  synthetic  noise  are  shown  in  Fig.  29.  Comparing  the 
actual  noise  in  Fig.  24  and  the  generated  noise  in  Fig.  29, 
the  generated  noise  correctly  models  the  random,  impulsive 
character  of  the  actual  noise  and  also  its  periodicity.  The 
synthetic  noise,  however,  does  not  incor-porate  the  constant 
bias  above  the  0  level  seen  in  the  actual  noise. 


F.   SUMMARY 

This  chapter  has  presented  a  detailed  analysis  of  the 
time-  and  frequency-domain  behavior  of  corona  noise.  A 
model  for  synthetic  corona  noise  was  postulated  based  on  a 
filtered  non-homogeneous  Poisson  process.   A  methodology  for 
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Figure  29.   Synthetic  Corona 

estimating  the  model  parameters  was  developed.  The  ad- 
vantage oF  this  model  over  previous  work  is  that  it  allows 
an  accurate  representation  of  the  time-domain  behavior  oF 
corona.  The  disadvantage  of  this  model  is  that  it  is 
scenario  dependent  and  -For  each  interference  situation  a  new 
set  of  model  parameters  must  be  specified.  One  limitation 
was  introduced  into  the  model  by  the  assumption  leading  to 
Eqn.  4-2.  The  lowpass  equivalent  -Filter  is  a  Function  oF 
the   bandpass    Filter   and    the   process    intensity,   and 
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consequent; ly ,  the   model  is  not  as  general  or  as  accurate  as 
the  gap  noise  model. 
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V .   ENERGY  DETECTION  RECEIVERS 

A.  INTRODUCTION 

Energy  detection  receivers  in  the  high  frequency  (HF) 
radio  band  operate  in  an  impulsive  radio  noise  environment. 
To  predict  and  improve  the  performance  of  energy  detection 
receivers,  accurate  models  of  the  signal  and  noise  envir- 
onment must  be  used.  The  noise  models  developed  in  earlier 
Chapters  will  be  used  to  simulate  atmospheric  and  man-made 
radio  noise  environments.  Two  types  of  energy  detection  re- 
ceivers will  be  considered:  a  -Past  Fourier  transPorm  (FFT) 
processing  system  and  a  compressive  receiver.  To  improve 
the  performance  o-P  these  receivers  the  concept  oP  the 
locally  optimum  receiver  is  introduced  to  suggest  a  robust 
post-detection  processor. 

B.  ENERGY  DETECTION 

The  energy  detector  receiver  measures  the  energy  in  a 
signal  over  a  specific  time  interval.  For  this  chapter  we 
will  consider  the  signal  to  be  a  modulated  sinusoid,  of  un- 
known frequency  at  the  receiver.  The  signal  is  then 
represented  as 

s(t)  =  Acos(u)0t  +  8(t))  (5-1) 

where  A  and  0)o  are  random  variables  and  9(t)  is  a  slowly 
varying  function  so  that  the  signal  power  is  effectively 
confined  to  a  narrow  bandwidth  B. 
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De-Pining  the  problem  as  a  hypo-thesis  test 

H0  :    x(t)  =  n(t)   noise  alone 

H, :  x(t)  =  n(t)  +  s(t)  signal  plus  noise.  (5-2) 
The  well  known  energy  detection  receiver  CRe-F.  42]  uses  te  o  , 
the  received  energy  as  the  test  statistic  such  that 

T 


Under  He 


n  (t)dt 


0 


Under  H,  :   tE 


Cn(t)  +  s(t) ]  dt 


(5-2) 


0 


For  a  sampled  bandpass  process  the  test  statistic  is 

N      2      2 
H0  :     tE  o  =  E  Cnc  .   +  n3  ,   ] 
1  =  1 


H, 


N 

=  E  Cnc  t 
i  =  l 


i  2    +    [n« 


,  : 


(  5  -  3  ) 


and 


xc  t      =  nc  t  +  Sci      or    xc  i  =  nc  i 
and  likewise  -For  the  quadrature  component. 

Fig.  30  shows  the  block  diagram  o-F  a  receiver  that  cal- 
culates this  statistic.  This  receiver  will  be  called  a 
square  and  sum  receiver.  The  square  and  sum  receiver  is  the 
optimum  receiver  -For  detecting  random  phase  signals  in 
Gaussian  noise  at  small  signal  levels.  At  large  signal  to 
noise  ratios  (SNR)  this  detector  is  very  close  to  optimum 
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Figure  30.   Square  and  Sum  Receivet 


CRef.  9:p.  3703.  When  N  =  BT  and  the  signal  is  sampled  at 
the  Nyquist  rate,  the  test  statistic  is  the  energy  in  the 
signal  in  a  period  T. 

When  it  is  desired  to  implement  the  receiver  shown  in 
Fig.  30  over  a  wide  frequency  band  relative  to  the  target 
signal  modulation  bandwidth  B,  a  number  of  options  may  be 
pursued.  The  channelized  receiver  approach  in  which  an 
independent  receiver  operates  at  each  -Frequency  will  often 
be  too  expensive.  Two  other  possible  methods  for  generating 
the  test  statistic  are;  (1)  a  compressive  receiver  as  a 
spectrum  analyzer  CRef.  43]  and  (2)  a  fast  Fourier  transform 
(FFT)  signal  processor. 
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It  has  been  shown  that  an  FFT  signal  processor  exactly 
implements  the  test  statistic  given  in  Eqn.  5-3  -For  the 
assumed  signal  in  a  single  FFT  filter  [Ref.  44].  Whether 
successive  samples  will  be  independent  or  not  depends  on  the 
number  of  trans-Forms  per  second  relative  to  the  signal 
bandwidth  B.  I-F  the  sample  interval  is  greater  than  1/B 
then  successive  samples  will  be  independent. 

In  the  compressive  receiver,  shown  in  Fig.  31,  a  -Fast 
sweeping  local  oscillator  linearly  scans  the  band  oF 
interest.  Any  narrow  bandwidth  signal  in  the  band  will 
appear  at  the  mixer  output  as  a  linear,  -Frequency  modulated 
(FM),  chirp  signal.  This  chirp  signal  is  passed  through  a 
weighting  Filter  to  reduce  sidelobes  and  into  a  dispersive 
delay  line  (DDL)  which  is  matched  to  the  inverse  oF  the 
linear  FM  sweep.  The  output  of  the  DDL  is  envelope  detected 
and,  when  referred  to  the  sweep  time,  it  provides  an 
estimate  of  the  spectrum  of  the  input  signal.  It  is  easily 
shown  that  for  a  sinusoid,  the  detector  output  is 
proportional  to  the  signal  magnitude.  Thus,  the  compressive 
receiver  approximates  the  square  root  of  the  ith  value  of 
the  test  statistic  in  Eqn.  5-3.  The  samples  will  be 
independent  if  the  scan  revisit  time  is  greater  than  1/B. 

C.   SQUARE  AND  SUM  RECEIVER 

In  this  Section  the  energy  detection  performance  of  the 
square  and   sum  receiver,  shown  in  Fig.  30,  will  be  examined 
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as  "the  noise  statistics  depart  -Prom  a  Gaussian  distribution. 
The  square  and  sum  receiver,  which  is  optimum  for  small 
signals  in  Gaussian  noise,  will  be  used  as  the  reference 
receiver  so  that  the  performance  of  improved  receiver 
structures  may  be  evaluated. 

The  number  of  samples,  N,  chosen  for  this  study  was  10. 
This  is  a  small  sample  size  and,  for  impulsive  noise,  the 
resulting  test  statistic  will  have  a  distinctly  non-Gaussian 
distribution.  In  Fig.  32  the  probability  of  false  alarm 
(PF» )  is  plotted  for  unit  root-mean-square  (RMS)  Gaussian 
noise  and  the  Hall  models  with  V0  equal  to  2.1,  3,  5  and  12. 
(see  Table  II).  The  normalized  threshold  is  in  reference  to 
the  unit  RMS  power  of  the  noise. 

A  Monte  Carlo  simulation  technique  was  used  to  obtain 
the  data  for  the  plot.  Since  the  data  aire  independent 
Bernoulli  trials  with  a  large  number  of  trials,  the  normal 
approximation  to  the  binomial  was  used  to  determine  the 
confidence  interval.  For  estimates  of  PF A  greater  than  5  x 
10"3  with  20,000  trials,  the  true  value  of  PF A  will  be 
within  20  percent  of  the  estimate  with  a  confidence  of  95 
percent  CRef.  55 : p .  282].  The  noise  samples  were  generated 
using  the  inverse  method  CRef.  45 :  p .  95.13  and  the  inverse 
cumulative  distribution  functions  (CDFs)  in  Table  I. 

In  Fig.  33  the  probability  of  detection  ( PD  )  curves  air& 
plotted  for  PF A  equal  to  1x10" a  .  A  rejection  method  CRef. 
45 : p .  9523   was  used   to  generate   the  random  deviates  since 
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the  CDFs  of  the  signal  plus  noise  densities  were  not  -Found. 
The  Pd  is  plotted  as  a  function  of  SNR  for  the  same 
parametric  conditions  as  used  in  Fig.  32.  As  anticipated 
the  Po  decreases  dramatically  as  the  noise  becomes  more 
impulsive.  In  the  case  where  V0  equals  12,  the  PD  decreases 
-Prom  greater  than  99  percent  to  1  percent  at  a  5  dB  SNR  for 
the  same  PF  a  . 

A  receiver  that  is  optimum  in  Gaussian  noise  suffers  a 
major  degradation  when  the  noise  statistics  become  non- 
Gaussian.  One  possible  approach  to  improving  the  per- 
formance is  to  design  a  receiver  that  is  based  on  the 
statistical  parameters  of  the  noise.  A  parametric  approach 
will  be  discussed  in  the  next  section. 

D.   LOCALLY  OPTIMUM  RECEIVER 

The  concept  of  a  locally  optimum  receiver  provides  a 
general  methodology  for  nonlinear  receiver  design  in  the 
presence  of  non-Gaussian  noise.  The  approach  is  to 
determine  a  receiver  structure  that  is  optimized  for  the 
difficult  small  signal  case  and  then  examine,  and  if 
necessary,  modify  the  structure  for  the  large  signal  case. 
In  this  Section  the  locally  optimum  receiver  for  energy 
detection  in  the  presence  of  Hall  type  impulsive  noise  will 
be  determined  under  the  assumption  of  weak  signal  levels. 
Previous  work  has  considered  known  signal  detection  in 
Cauchy  noise   which  is   related  to   one  of   the  Hall   models 
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[ReF.  46].  AFter  optimizing  the  design  -For  weak  signal 
levels  the  performance  oF  the  receiver  will  be  examined  -For 
moderate  to  large  signal  levels.  The  receiver  structure 
will  be  parametric;  that  is  it  will  depend  on  the  noise 
density  -Function  parameters. 

To  illustrate  the  method  For  deriving  the  locally 
optimum  receiver,  consider  the  Following  binary  hypotheses 
test.  We  are  given  a  sequence,  xt  ,  o-P  N  samples.  Under  the 
null  hypotheses  H0  ,  we  assume  that  the  sequence  consists 
only  of  independent  and  identically  distributed  noise 
observations  with  a  common  density  p(n) .  Under  the 
alternative  hypothesis  H, ,  we  assume  that  the  sequence 
consists  oF  an  additive  mixture  o-F  a  small  signal  and  noise. 
The  problem  may  be  stated  as 
Ha  :    x,  =  n,  noise  alone 

Hi :    x,  =  n(  +  As,         signal  plus  noise   1  =  1 , . . ,N 
with  A  an  unknown  small  amplitude. 

The  optimum  receiver,  in  terms  oF  maximizing  the 
probability  oF  detection  while  keeping  the  probability  oF 
False  alarm  below  a  certain  level  is  given  by  the  Neyman- 
Pearson  lemma.  This  receiver  compares  the  likelihood  ratio 
test  statistic  to  a  threshold  which  is  chosen  to  achieve  the 
desired  probability  oF  False  alarm. 

The  liklihood  ratio  test  is 


A(x;A)  = 


Po  <x|H0  ) 


>  T  decide  Hc 


<  T  decide  Ht 


or 


<:  5  -  4  ) 
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where  p0 (x|HQ )   and  pi (x|H, )   are   the   probability   density 
-Functions  for   the  observation   vector  x   under   Ha   and   H, 
respectively.    Since  the   noise  samples  ar&     independent  and 
i  dent ically   di  str i  buted ,    pQ  ( x  |  Ha  )   and   pt  (x | Ht  )   can   be 
written  as  products  of  the  univariate  noise  densities: 


and 


N 

p«  (.*.)  =  n  pt  (x, ) 

i  =  l 


N 
p ,  ( x )  =  n  p ,  (  x  t   -  As(  )  . 

1  =  1 


(5-5) 


(5-6) 


The  test  can  then  be  written 


N 


A  (  x  ;  A )  = 


n  pi  ( xt  -  As,  ) 
1  =  1 

N 

n  pt  (xt  ) 

i  =  l 


>  T  decide  HQ  ,  or 


<  T  decide  H, 


(5-7) 


In  deriving  the  locally  optimum  receiver  several 
approaches  have  been  used.  Rudnick  CRef .  473  derived  a  form 
of  the  locally  optimum  receiver  by  expanding  the  likelihood 
ratio  in  a  Taylor  series  about  x  and  truncating  after  the 
second  order  terms.  Capon  CRef.  48]  showed  that  the  locally 
optimum  detector  maximizes  the  slope  of  probability  of 
detection  versus  signal  strength  function  at  the  origin  and 
employs  a  statistic  based  on  the  following  test 

■oA  (  x  ;  A  ) 


tL  o  ( x )  = 


-oA 


(5-8) 


A  =>  0 
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In  Appendix   H  it   is  shown   that   the   locally   optimum 

receiver  -For  the  detection  o-f  a  random  sinusoid  is 

N 
t  \0   (r)  =    E  g(r,  )  (5-9) 

1  =  1 

where  g(- )   is  the   zero-memory  nonlinearity   (ZMNL)  defined 

be  low 

F'  '  (r)      f  (r) 

g(r)  =  +  (5-10) 

•P(r)       r-P(r) 

and  r   is  a  vector  oF  envelope  samples.   The  -Function,  -P  ( r )  , 

For  the  Hall  model  oF  atmospheric  noise  discussed  in  Chapter 

II  is 

9-1 
•  KY 

f  (  r )  =  (5-11) 

2     2  (8-D/2 
(r   +  y  > 

where  K   is  a  constant.   From  Eqn.  5-10  the  ZMNL  oF  the  test 

statistic  is 

2 
(9+1 ) (8+3)r         2(8+1 ) 

g  ( r )  =  -    .        (5-12) 

2    2  2  2     2 

( r  +  y  >  ( r   +  y  > 

The  receiver  that  implements  this  test  is  shown  in  block 
diagram  Form  in  Fig.  34.  The  ZMNL's  For  the  unit  root-mean- 
square  Hall  models  (see  Table  II)  parameterized  by  VD  are 
shown  in  Figs.  35  and  36.  The  ZMNL  receiver  characteristic 
For  values  oF  VD  equal  2.1,  3,  6,  10  and  14  are    plotted. 

Recalling  that  the  nonlinearity  is  optimized  For  van- 
ishingly  small  signals  the  ZMNL  characteristic  may  be 


109 


<rt  > 


Yes 


No 


H0 
Figure  34.   Locally  Optimum  Receiver  for  Energy  Detection 

intuitively  analyzed.    Small  values  of  the  sampled  envelope 

will  be   given  a  negative  weighting,  intermediate  values  are 

emphasized  and   large  values   of  the   sampled   envelope  are 

effectively  nulled.    As   the  noise   becomes  more  impulsive, 

which  corresponds  to  a  larger  value  of  VB  and  smaller  9,  the 

above  features  become  more  pronounced.   This  is  seen  in  Fig. 

35  where   the  ZMNL  has  large  negative  values  that  extend  off 

the  plot   for  very  small  envelope  values.   The  maximum  value 

of  the  ZMNL  occurs  at 

1/2 
x    =    y( (Q+5) / (Q+l) )  (5-13) 

and  the  zero  crossing  at 

1/3 
x  =  y(2/ (9+1 ) )     .  (5- 1 4) 

Both  of  the  points  are    proportional  to  the  mean  value  of  the 

underlying  distribution  for  the  same  Hall  model  number. 
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Fig.  37  is  the  PFA  versus  threshold  for  the  locally 
optimum  receiver  shown  in  Fig.  34.  The  number  N  of  samples 
is  equal  to  10.  The  Hall  noise  samples  are  parameterized 
with  the  same  VD  values  used  in  Fig.  32.  The  threshold  data 
obtained  -From  this  plot  was  used  to  generate  the  P0  curves 
shown  in  Fig.  38.  The  only  complete  data  set  was  obtained 
For  the  Hal  12  case  with  V0  equal  to  2.1.  Some  data  was 
obtained  For  V&  equal  to  5.  The  peaks  in  P0  For  the  more 
impulsive  noise  occur  at  smaller  SNR  levels,  not  shown  on 
the  plot.  The  data  points  are  plotted  on  the  graph  since 
the  dashed  curves  only  approximate  the  true  values. 

It  can  be  seen  that  the  locally  optimum  receiver  does 
improve  the  Po  (For  the  same  PF a )  in  the  small  signal  case. 
In  comparing  the  values  with  Fig.  33  For  V0  equal  to  2.1  the 
locally  optimum  receiver  has  a  performance  better  than  the 
square  and  sum  receiver  From  -10  dB  to  7  dB  SNR.  Above  7  dB 
SNR,  the  Po  For  the  locally  optimum  receiver  drops  to  0. 
This  is  due  to  the  Fact  that  the  receiver  is  optimized  For 
small  signals.  The  Hall2  ZMNL  in  Fig.  36  shows  the  reason 
For  this.  It  peaks  at  approximately  1.3  times  the  RMS  value 
o-P  the  noise  and  approaches  0  as  the  input  goes  to  infinity. 
This  behavior  of  the  ZMNL  must  be  modified  to  obtain  a 
receiver  that  will  also  work  in  the  large  signal  case.  In 
the  next  Section  a  robust  receiver  structure  is  proposed 
that  does  this.  A  robust  receiver  is  differentiated  From  a 
parametric  receiver   in  that   the  robust   receiver  does   not 
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Figure  37.   PF  A  -For  Locally  Optimum  Receiver  in  Hall  Noise 
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depend  on  the  parameters  of  the  noise  density  and  it  per- 
forms well,  i-P  subopt  1  ma  1  1  y  ,  over  a  class  or  even  differing 
classes  of  noise  densities. 

E.   ADAPTIVE  LIMITER 

The  ZMNL  derived  in  the  previous  Section  is  a  relatively 
complex  function.  Furthermore,  it  does  not  work  for  large 
signal  levels.  In  this  section  we  consider  a  practical 
implementation  of  the  ZMNL  consisting  of  an  adaptive  limiter 
to  approximate  the  ZMNL  function  with  a  modification  to 
improve  large  signal  performance.  Although  the  locally 
optimum  receiver  was  derived  for  Hall  type  impulsive  noise, 
if  the  approximation  to  it  is  properly  designed,  it  should 
show  robust  performance  in  the  presence  of .many  impulsive 
noise  process. 

Examining  the  Po  curves  of  the  locally  optimum  receiver 
for  the  Hal  12  model  shows  that  it  provides  increased 
detection  performance  over  the  square  and  sum  receiver  from 
-10  dB  to  7  dB  SNR.  Since  this  region  is  where  the  square 
and  sum  receiver  suffered  the  worst  performance  degradation 
in  impulsive  noise,  the  Hall2  ZMNL  was  chosen  to  be  modified 
to  obtain  a  robust  non 1 i near i ty .  The  solid  curve  in  Fig.  39 
is  the  Hal  12  ZMNL.  To  correct  its  performance  at  large 
signal  levels,  the  ZMNL  was  first  altered,  as  shown  by  the 
dashed  line  in  Fig.  39,  by  extending  its  characteristic 
straight  out   from   its   maximum   value.      Furthermore,   an 
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equivalent  test  statistic  is  obtained  if  g'Cx)  is  any  linear 
function  of  g(x)  CRef.  463  so 

g'Cx)  =  ag(x)  +  0    a     \    <2 .  (5-15) 
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Figure  39.   Modification  to  ZMNL 


With  suitable  parameters  a     and  fl,  the  modified  ZMNL  is  shown 
in  Fig.  40.   The  level,  a,  remains  to  be  determined. 

Recalling  Eqn.  5-13,  the  argument  of  the  ZMNL  char- 
acteristic at  its  maximum  value  is  proportional  to  the  mean 
value  of  the  Hall  distribution  for  which  it  was  derived  (see 
Table  I).  For  the  Hall2  distribution  the  proportionality 
constant  is  approximately  1.6.  Therefore,  in  order  to 
preserve  the  improvement  in  P0  in  the  same  SNR  region  as  the 
Hall2  ZMNL,  the  limiter  level  was  set  to  equal  1.6  times  the 
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g(x) 


Figure  40.   Linearized  Modification  to  ZMNL 

mean  value  o-P  the  assumed  noise  distribution. 

A  block  diagram  o-P  the  receiver  is  shown  in  Fig.  41 
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Figure  41.   Receiver  Structure  -For  Adaptive  Limiter 
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Note  that,  the  ZMNL  must  be  adjusted.  This  information  could 
be  known  a  priori  or,  as  symbolically  shown  in  Fig.  41, 
determined  -Prom  the  statistics  o-P  the  noise.  Both  cases 
will  be  examined  below. 

The  -False  alarm  rate  performance  o-P  this  modified 
receiver  in  Hall  and  Gaussian  noise,  where  the  statistics 
are  known  a  priori,  is  shown  in  Fig.  42  Por  N  equal  to  10. 
The  value  Por  a  was  determined  Prom  the  noise  mean.  At  a 
PF  a  equal  to  1  x  10"  a  ,  the  P&  curves  o-P  this  receiver  are 
shown  in  Fig.  43.  The  threshold  was  determined  -Prom  Fig. 
42.  The  Gaussian  noise  perPormance  o-P  this  receiver  is 
poorer  than  that  o-P  the  square  and  sum  receiver  but  only  by 
a  -Pew  percent.  However,  there  is  a  dramatic  improvement  in 
the  Pi>  For  the  impulsive  noise  cases.  It  is  interesting  to 
note  that  Por  the  same  noise  power,  a  well  designed 
nonlinear  receiver  provides  better  detection  performance  in 
impulsive  noise  than  the  optimum  receiver  for  Gaussian  noise 
CRef.  49]. 

By  plotting  the  empirically  determined  threshold  from 
Fig.  42  versus  the  mean  of  the  noise  distribution  it  is  seen 
that  the  relationship  is  almost  linear.  This  is  shown  in 
Fig.  44  for  two  false  alarm  rates.  Using  this  relationship 
the  ZMNL  and  threshold  can  be  adaptively  modified  to  provide 
nearly  constant  false  alarm  rate  (CFAR)  performance  under 
varying  noise  conditions.  This  is  a  desirable  charac- 
teristic for  an  energy  detection  receiver. 
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Figure  42.   PP A  For  Limiter  Receiver  in  Gauss  and  Hall  Noise 
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A  -Pinal  question  to  be  resolved  in  implementinq  the 
receiver  is  the  real  time  determination  o-P  an  estimate  01-' 
the  mean  and  mean  square  values  o-P  the  noise.  One  approach 
is  to  sample  the  noise  in  a  signal  -Pree  adjacent  channel  and 
another  is  to  Porm  the  sample  mean  o-f  the  noise  using  past 
sample  values  with  no  signal  present  CReP.  50].  We  will 
implement  the  latter  approach  and  test  the  adaptive  receiver 
on  actual  signal  and  noise  samples.  In  sequencing  the 
samples  two  general  approaches  may  be  used.  The  samples  may 
be  processed  in  blocks  with  no  overlap  which  is  the  case 
corresponding  to  the  hypothesis  tests  discussed  in  this 
Chapter.  The  other  option  is  sequentially  processing  each 
sample  with  rank  two  updates  oP  the  test  statistic,  the  ZMNL 
and  the  threshold.  In  this  case  the  hypothesis  test 
Pormulation  will  provide  a  lower  bound  on  P&  . 

The  sequential  sampling  approach  was  chosen  and  two  runs 
oP  the  adaptive  receiver  ar&  shown  in  Figs.  45  and  46.  The 
test  statistic,  shown  by  the  dotted  line,  was  Pormed  Prom 
the  10  previous  samples  using  the  receiver  structure  in  Fig. 
41.  The  adaptive  threshold,  shown  by  the  dashed  line,  was 
Pormed  Prom  the  sample  mean  and  mean  square  estimates  oP  the 
previous  20  though  50  samples.  The  sample  mean  is  shown  by 
the  heavy  solid  line  and  the  individual  envelope  data  points 
by  the  light  solid  line. 

This  receiver  structure  requires  H0  to  be  true  at  the 
start  oP   the  test   and  will   detect  the   transition  to   H, . 
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Since  the  adaptive  threshold  assumes  noise  only  samples,  it 
will  be  corrupted  when  signal  plus  noise  samples  ar&  used  to 
■Form  it.   This  will  occur  20  samples  after  Hi   becomes  true. 

This  effect  can  be  seen  in  Fig.  45  where  an  actual 
signal  -From  one  Filter  oF  a  compressive  receiver  is  plotted. 
The  noise  is  highly  impulsive  and  would  be  a  very  difficult 
detection  environment  For  a  receiver  optimized  For  Gaussian 
noise.  Samples  50  to  110  are  noise  only  and  samples  110  to 
220  are  signal  plus  noise.  The  test  statistic  exceeds  the 
threshold  at  sample  125  and  the  receiver  would  declare  H, 
true  at  this  point.  However,  due  to  the  Fact  that  the  test 
statistic  and  the  threshold  are  both  corrupted  by  the  signal 
plus  noise  samples,  the  test  statistic  later  drops  below  the 
new  threshold. 

In  Fig.  46,  the  receiver  performance  in  a  simulated  gap 
noise  environment  is  analyzed.  The  signal  commenced  at  sam- 
ple 110  and  the  receiver  declared  the  signal  present  at 
sample  117.  Again  the  threshold  is  later  corrupted  and  the 
test  statistic  drops  below  the  threshold. 

F.   SUMMARY 

This  Chapter  has  presented  a  systematic  development  oF 
the  bandpass  energy  detection  problem  in  impulsive  noise. 
Starting  with  a  Hall  density  Function  For  the  noise,  the 
locally  optimum  receiver  is  derived.  This  receiver  is 
modified   to    design   a    practical   receiver   with   robust 
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performance  in  non-Gaussian  noise.  The  design  methodology 
clearly  delineates  the  role  of  the  test  statistic,  threshold 
and  nonl mear ity .  This  allows  independent  evaluation  of 
each  of  these  features.  One  aspect  of  energy  detection 
receivers  not  considered  is  the  role  of  data  in  filters 
adjacent  to  the  filter  being  processed.  In  an  impulsive 
noise  environment,  the  noise  statistics  of  adjacent  filters 
will  be  correlated  and  their  use  may  potentially  offer  a 
further  performance  gain. 
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V I •  CONCLUSION 
A.   RESULTS 

The  dissertation  has  investigated  the  modeling  oP 
atmospheric  and  man-made  radio  noise  and  applied  the  results 
to  the  bandpass  energy  detection  problem. 

A  generalized  means  oP  describing  an  impulsive  noise 
process  in  terms  o-P  its  complex  envelope  was  derived. 
Previous  models  oP  atmospheric  noise  were  surveyed  and  some 
extensions  to  the  Hall  model  Por  atmospheric  noise  were 
developed . 

Using  the  physical  characteristics  oP  gap  noise,  a  nine 
parameter  model  was  constructed  that  allows  arbitrary  sour- 
ces oP  gap  noise  to  be  synthesized.  The  driving  impulses 
Por  the  gap  noise  were  postulated  to  be  points  oP  a 
branching  renewal  process.  The  application  oP  this  type  o-P 
probabilistic  model  to  power  line  noise  is  unique. 
Additional  results  include  a  new  prooP  oP  the  power  spectrum 
oP  an  equilibrium  renewal  point  process  and  a  derivation  oP 
the  power  spectrum  oP  a  truncated  branching  renewal  point 
process.  The  postulated  model  was  shown  to  correctly 
predict  all  the  significant  Peatures  oP  the  observed 
spectrum . 

In  a  similar  Pashion,  corona  noise  was  postulated  as  a 
Pi ltered  impulse  process  with  a  periodically  modulated  rate 
Punction.    This  model   was  shown  to  correctly  duplicate  the 
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observed  "time  domain  behavior  of  corona  noise.  A  method  for 
estimating  the  parameters  of  the  model  was  discussed  and  an 
actual  example  of  corona  noise  simulated. 

The  energy  detection  problem  was  also  considered  and 
starting  with  an  assumed  density  for  impulsive  noise,  a 
design  methodology  for  a  parametric  and  robust  receiver  was 
developed.  This  receiver  was  then  tested  against  actual  and 
simulated  data  and  shown  to  be  superior  to  the  Gaussian 
noise  receiver  when  corrupted  with  impulsive  atmospheric  or 
power  line  noise. 

It  has  been  stated  that  in  engineering  design  one  seeks 
not  so  much  to  be  optimum  but  to  avoid  crippling  non- 
optimalities  CRef.  513.  It  is  hoped  that  this  research  will 
allow  system  designers  to  test  their  systems  with  simulated 
interference  based  on  the  models  presented.  The  robustness 
of  the  system  in  an  actual  noise  environment  can  then  be 
evaluated,  as  was  done  in  this  research  for  the  locally 
optimum  energy  detector.  The  noise  models  developed  here 
should  be  particularly  applicable  to  systems  that  operate 
from  a  fixed  site  within  line  of  sight  of  power  lines  where 
the  chance  of  having  interference  from  gap  or  corona  noise 
sources  is  significant. 

B.   FURTHER  RESEARCH 

At  the  conclusion  of  this  dissertation  several  problems 
ar&    worthy   of  further   study.    Research  should   be  done  to 
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de-Pine  the  range  of  values  for  the  power  line  noise  model 
parameters  and  to  determine  the  effect  of  power  line 
construction  practices  on  the  parameters.  The  suscept- 
ibility of  various  types  of  communication  equipment  to 
simulated  gap  and  corona  noise  should  also  be  examined.  The 
corona  noise  model  needs  to  be  extended  so  that  it  can  be 
specified  and  estimated  using  a  bandpass  vice  lowpass  filter 
impulse  response.  The  work  on  energy  detection  suggests 
many  opportunities  for  further  research.  The  problem  of 
estimating  the  unknown  noise  parameters  for  the  adaptive 
updates  and  the  inclusion  of  adjacent  filter  samples  into  an 
algorithm  aire    prime  examples. 
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APPENDIX  A 
INSTRUMENTATION 

The  instrumentation  configuration  employed  to  provide 
data  on  the  detailed  time-  and  -Frequency-domain  properties 
oF  high  -Frequency  (HF)  radio  noise  is  shown  in  Fig.  47.  The 
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Figure  47.   Simpli-Fied  Block  Diagram  o-F  Measurement  System 

radio  Frequency  (RF)  input  to  the  instrumentation  was  either 
a  Fixed  HF  long  wire  antenna  or  a  whip  antenna  mounted  on  a 
mobile  van.  One  complete  measurement  system  is  trans- 
portable and  may  be  moved  From  site  to  site.  The  other 
system  is  installed  in  a  mobile  van  For  noise  measurements 
directly  at  the  interference  source  CReF.  53]. 

In  the  lower  halF  oF  Fig.  47,  a  Hewlett-Packard  141T 
Spectrum  Analyzer  is  used  as  a  scanning  receiver  to  drive  a 
Develco  7200B  3-  Axis  Display.   The  spectrum  analyzer  can  be 
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■tuned  to  any  desired  frequency  in  the  HF  band.  Its  scan 
rate,  scan  width,  intermediate  frequency  (IF)  bandwidth,  IF 
gain,  RF  attenuation,  and  other  controls  can  be  adjusted  to 
best  describe  the  noise  under  observation. 

An  alternate  and  complementary  part  of  the  system  is 
used  to  examine  the  narrowband  properties  of  HF  noise  and  is 
depicted  in  the  upper  half  of  Fig.  47.  An  HF  receiver  is 
used  as  an  amplifier/translator  and  tuned  to  a  frequency 
where  the  noise  is  present.  The  demodulated  audio  output  is 
then  applied  to  a  Wavetek  (Nicolet)  UA500A  spectrum  analyzer 
which  subdivides  the  audio-output  spectrum  into  about  500 
frequency  segments.  When  the  full  audio-output  bandwidth  of 
5  kHz  is  examined,  the  UA500A  provides  a  frequency 
resolution  of  10  Hz  per  segment.  The  Wavetek  analyzer  can 
provide  individual  transforms  or  average  a  selected  number 
of  successive  transforms.  The  transforms  (either  individual 
or  averaged)  are  then  presented  on  the  3-axis  display  or  ar\ 
osci 1 loscope. 

The  3-axis  display  provides  a  continuous  moving  real- 
time visual  representation  of  the  analog  output  from  the 
spectrum  analyzers  or  narrowband  receiver  audio  output.  The 
receiver  or  analyzer  output  analog  data  are  digitized  in  the 
display  and  stored  in  a  semiconductor  memory.  The  data  in 
the  memory  3ir&  formatted  and  shown  on  a  cathode  ray  tube 
(CRT)    in   a   convenient   f requency-ampl l tude-t l me   (3-axis) 
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■Format.  The  3-axis  presentation  can  be  -Frozen  at  any  desired 
time  and  photographed  with  a  standard  oscilloscope  camera. 

Fig.  48   illustrates  the   procedure  used   to  make  the  3- 
axis  presentation.   The  analog  input  is  divided  into  512 
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Figure  48.   Diagram  o-P  Data  Format 

equally  spaced  data  points  (indicated  by  the  horizontal  dots 
oF  Fig.  48).  The  signal  amplitude  at  each  data  point  is 
represented  by  an  8-bit  word.  When  a  scan  is  completed,  its 
data  are  also  stored  in  memory.  Line  1  in  the  view  moves  to 
line  2,  and  the  new  scan  appears  as  line  1.  Subsequent 
scans  move  earlier  data,  line  by  line,  upwards  along  the 
time  axis  to  create  a  rising  raster  type  display.  When  the 
memory  is  Full  (64  scans),  each  new  scan  is  entered  into  the 
bottom  line,  and  the  oldest  data  at  the  top  line  is 
discarded.  The  resulting  animated  view  provides  a  visual 
picture  oF  noise  and  signals  within  the  block  o-F  Frequencies 
lF  the   spectrum  analyzer   output  is  being  observed.    IF  the 
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receiver  audio  output  is  being  observed  the  display  is  a 
stacked  series  o-F  consecutive  time  records. 

The  3-axis  display  has  a  number  oF  controls  to  aid  the 
observer  in  analyzing  the  structure  oF  data  presented. 
Among  these  controls  are:  (1)  elevation  and  azimuth  geometry 
controls  to  vary  the  viewing  aspect,  (2)  an  amplitude 
compression  control,  (3)  a  threshold  control  to  \/ary 
background  levels,  (4)  time-axis  expansion  controls,  and  (5) 
a  stop-action  switch  to  -Freeze  the  data  in  memory  -For 
detailed  observation,  For  observation  From  various  aspects, 
or  For  photographing.  These  controls  are  used  to  optimize 
the  presentation  o-F  desired  signal  detail. 

Accurate  Frequency,  time  and  amplitude  calibrations  ar& 
maintained  so  that  the  resulting  3-axis  photos  can  be 
manually  scaled  For  precise  signal  detail.  The  digitized 
data  at  the  display  input  ar-&  available  at  a  digital 
interface  connector  For  external  digital  recording  or 
processing.  Received  data  can  be  processed  by  computer  to 
obtain  conventional  noise  statistics.  This  Feature  permits 
the  comparison  oF  selected  3-axis  views  oF  noise  with 
standard  statistical  descriptors  oF  noise. 
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APPENDIX  B 
FILTERED  IMPULSE  PROCESS 

The  complex  envelope  of  a  noise  process  can  be  described 

by  the  following  equation 

N(t)   -j8i 
n(t)  =  g(t)  +  E  a, e     h(t-tt)  (B-l) 

i  =  l 

where  g(t)   is  a   quadrature  Gaussian   process,   h(t)    is   a 

linear  -Filter,  at   is  an  independent,  identically  distributed 

random   amplitude    and   9,    is   a   random   phase   uniformly 

distributed  over  0  to  2n  .   The  joint  characteristic  -Function 

of  the   inphase   and   quadrature   components   of   the   above 

process  is 

j  (u)t  nc  (t)  +  oia  na  (t)  ) 
0  ((Ut  ,o)a  )  =  ECe  ]  .         CB-2) 

IF   the    Gaussian   noise   is   assumed   independent   oF   the 

impulsive  noise   the  characteristic  Function  of  the  Gaussian 

term  alone  can  be  written 

j  (ait  gc  (t)   +  u)a  g»  Ct)  ) 
0a  (on  ,(ja  )  =  ECe  ]  .      (B-3) 

The  inphase   and  quadrature   components   oF   bandpass   white 

Gaussian  noise  a\r&     independent  so 

2      2 
-NoBCWt    +  Wa   )/2 
09  (Ut  ,o)a  )=  e  .  (B-4) 

The  joint   characteristic  Function   oF  the   ith  pulse  oF  the 

inphase  and  quadrature  impulsive  noise  terms  is 
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jai   ((Ui   (hc  (t-t,  )cos8i      -    h8  (t-tk  )sin8i  )     + 
0  t    (  0Jt    ,  Ida  )     =    E  C  e 

(i)a  (hc  (t-t»  )sin9i      +    h3  (t-tt  )cos6,   )) 

].  (B-5) 

Since  the   noise  process   is  narrowband,   its  joint   density 

■Function   is    circularly   symmetric    and    the    -Following 

transFormat ion  can    be  applied 


2      2  1/2 


-1 


a  =  tan   ( ojt  /oja  ) 
u)i  =  o)r  cosoc  cja  =  (DrSina.         (B-6) 

This  type  oF  transFormat ion  oF  a  2-D  Fourier  transForm  to  a 
one  dimensional  transForm  oF  the  envelope  is  known  as  a 
Hankel  TransForm  and  the  inversion  Formula  is 


p(r)  = 


rajr0((j^  )  J0  (ruJr  )  du)r 


(B-7) 


0 


where  J0    is  the  ordinary  Bessel  Function  oF  the  First  kind, 
order  0.    To   simplify  the   derivation  it   is  assumed   that 
h9(t)   is    zero   (see    comment   before    Eqn.   2-21).   Now 
transForming  Eqn.   B-5  and   taking  the   expected  value  oF  8, 
and  a, 

co   2n 
0,  (air  )  =  EC 

0   0  (B-8) 

Since  9,   is  uniformly  distributed  over  0  to  2n , 


jairai  hc  (t-ti  )cos(9t  -a) 
p (9, )p(a, ) e  d9  ,  da,  ] 
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0t  (ti>P  )  =  EC   p  Cat  )  Jo  (uir  a i  hc  Ct  )  )  dat  3  . 

0  (B-9) 

If  N(T)  in  Eqn.  B-l  is  a  Poisson  process  with  rate  X  and 
observation  interval  T,  the  t,  's  are  independent  and  uniform 
over  (0,  T)  allowing  the  characteristic  -function  of  the 
impulsive  term  to  be  written  as 


<d   k 
0i  (ur  )  =   E  0  (ajr  )  PrCN  =  kD  . 
k=0 


CB-10) 


Now  using   the  probability   law  for   k  events   of  a   Poisson 
process 


a>   k       -XT     k 
0,  <u>,.)  =  Z    0  Cu)r  )  e    (XT)  /k  ! 
k=0 


(B-l 1 ) 


which  sums  to 


X 


p(a) 


[Jo  (ojrahc  (t))  -  lDdTda 


(B-12) 


0i  (o)^  )  =  e 

Using   the    convolution  property    of  zero   order   Hankel 

transforms  CRef.   54]  and  Eqn.  B-7   the  density  function  of 

the  envelope  for  Gaussian  plus  impulsive  noise  is 


p(r)  = 


rcjr  Jo  (ru)P  )0a  (<ur  )0i  ( cjr  )  du)r  . 


(B-13) 
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APPENDIX  C 

SIGNAL  PLUS  NOISE  DENSITY 

For  the   truncated  Hall   model  with   8  =   3,  the  joint 

density  -function   o-P  the   m-phase  and  quadrature  components 

is  given  by: 

2    2 
D   2y  2    2      2 

p(x,y)  =  0  <  x  +  y   <  TP 

2      2    2    2  2 
2n(D  -  1  )  <x  +  y  +  y  > 

2    2      2 
=   0  x+y>TP  CC-1) 

where 

2   2  1/2 
D  =  (  1  +  T„  /y  ) 

The  sum   of  a   random  phase  sinusoidal  signal  oP  amplitude  A 

and  noise   can  be   expressed  in   terms  o-P   its  ln-phase   and 

quadrature  components  w  and  z  where: 

w  =  AcosCiy)  +  x 

z  =  Asin(iy)  +  y.  (C-2) 

Eqn.  C-2  is  substituted  into  Eqn.  C-l  and  transformed  to 
polar  coordinates  with 

w  =  rcos(0) 

z  =  rsin(0)  (C-3) 

where  r  represents  the  envelope  and  0  the  phase  o-P  the 
signal  plus  noise.  The  joint  probability  density  function 
p(r,0)  conditioned  on  uy  is 
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p  (r  ,0  |  vy)  = 


D  r2y 

2      2 
2n(D  -l)(r  -  2Arcos(i+/  -  0) 


2    2  2 
A  +  y  ) 


2  2      2 

0  <  r  -  2Arcos(\|/  -  0)  +  A   <  TP 

=  0  elsewhere.  (C-4) 

This  expression  is  valid  -For  all  values  of  r  between  0  and 
Tp  -  A.  Since  our  principal  interest  is  the  small  signal 
case,  A  <<  TP  ,  this  will  not  affect  subsequent  analysis.  By 
integrating  Eqn.  C-4  with  respect  to  0,  the  conditional  mar- 
ginal density  -Function  of  the  envelope,  p(rfi|/),  is  obtained; 


p  (  r  |  iy)  =  2 


p(r  ,0  |  \y)  d0  . 


(C-5) 


0 


The  integral  is  a  Function  only  oF  cos(iy  -  0 )  ,  \y  appears 
nowhere  else  in  the  integral.  Therefore  by  setting  a  =  \y  - 
0  and  using  Eqn.  2.554-3  in  [Ref.  543  gives 


2  2  2  2       2 

(r    +    A    +    y     >     2y    rD 


p  (r  |  4/)     = 


4  22  22  22442 

(r-2Ar       +    2A    y       +    2y    r    +    A    +    y     )(D    -     1) 


(C-6) 
Since  \y  is  assumed  independent  From  sample  to  sample  and 
uniformly  distributed  over  2n , 

2n 


p(r)  = 


p(r|  iy)p(vy)  d\|/ 


(C-7) 


0 


and 
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2  2  2  2       2 

Cr    +.  A    +    y     >     2y    rO 


4  22  22  22  442 

(r-2Ar       +2Ay       +2y»~2+A    +    Y><D-l) 

(C-8) 

By  a   similar  argument,   the  density   function  -For   a  random 

phase  sinusoidal   signal  in   narrowband  noise   with  a   Hal  15 

distribution  is 

4    2      2 
2y  r(b   +  2a  ) 

p(r)  =  — — —  (C-9) 

2     2  5/2 
(a   -  b  ) 

where 

2     2      2 

a    =    r       +  A   +   y 

and 

b  =  -2Ar. 


140 


APPENDIX  D 
SPECTRUM  OF  A  RENEWAL  PROCESS 

The  complex   envelope  representation  o-f  f  1  ltered  impulse 

noise  on  a  -Finite  interval  T  is  given  by 

N(T)   -j9t 
n(t)  =    Z    at  e  h(t-t,  )  ;         0  <  t  <  T. 

1  =  1 

CD-I  ) 

N(T)  is   a  random   unit  counting  process  denoting  the  number 

o-f  impulses  in  the  interval,  the  -Ctt  >  are  the  random  arrival 

times  and   8t   is   uniformly  distributed   over  0   to  2n .   The 

pulse    amplitudes     -Cat  >    ar&  identically    distributed 

statistically   independent    random   variables   with   second 

moment  Aa  and  -fourth  moment  A4  .    The  -filtering  e-f-fect  on  the 

impulses  ar&       accounted  -for   by  the   time  invariant   complex 

envelope  impulse  response  h(t). 

The  envelope  squared  o-f  n(t)  is 

2  N(T)  2 

E  (t)  =  n(t)n- (t)  =   E  at  h ( t-tt  ) h" ( t-tt  ) 

1  =  1 

(D-2) 

where  we   have  assumed  terms  o-f  the  -form  h(t-tt  )  h"  (t-tk  )  -for 

different  l   and  k  are    negligible.    We   estimate  the  power 

spectral   density    (PSD)   o-f    the   envelope   by   averaging 

magnitude  squared,   length  T,   Fourier  trans-forms   o-f  Ea  (t). 

The  mean  value  o-f  the  estimate  is 

S(u),T)  =  ECEa  (cj)Sa  -  (cu)  ]/T  (D-3) 
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where 


£a  <0),T)  = 


2     -jut     N(T)   2  -jart, 
E  (t)e     dt  =  E   a,  e      Ha  <o>) 
1  =  1 


0 


and 


Ha  (a>)  = 


-.idlt 


h(t)h-  (-t)e 


d-t 


is  -the   Fourier  transform   oF  the   magnitude  squared  impulse 
response.    We  have   assumed  the   impulse  response   is   much 
shorter  than  the  observation  interval  T. 
The  mean  oF  the  estimate  becomes 

_  N(T)  N(T)   2   2   -ju)(tt  -t„  )  2 

S(u),T)  =  EC  E     Z       at     a*        e  ]  |  Ha  (  id  )  |  /T. 

k=l   i=l 

(D-4) 

This  is   a  compact   and  general   expression  -For  the  PSD  o-F  a 

truncated  Filtered   point  process.     It   is  in   terms   oF   a 

summation    oF    the    characteristic    Functions    oF    the 

interarrival  time   distributions  between  all  combinations  oF 

pairs  oF   points.    (The   term   incorporating   the   Frequency 

response  oF   the  Filter   will  be   set  to  one  For  the  rest  oF 

the  analysis.)   Conditioning  on  N(T) 

«      N    N     2   2  -jti)(tl  -t„  ) 
S(id,T)  =  E   EC  E    E   a,  a„  e  ]PrCN  =  1]/T. 

1=2    k=l  i=l 

(D-5) 

Plotting  the   difference  k-j  indicates  a  different  summation 

over  m  and  n  shown  below. 
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Let    m       =     i        -     k       and    n    =     i     •+■    k    such    that    t.     =    tt      -    tk      is    the 
rnth     interarrival    time.        Then    Eqn.     D-5    can    be    written    as 


_  oo  N-l-n  2       -joit. 

S(d),T)     =       E       E-C    E  E       Aa     e  +    NA4 

1=2  n=l     m=l 


N- 1        n       2       -  jut. 

E  E    Aa     e  >PrCN    =     IJ/T. 

n=l     m=l 


(D-6) 


On  -Further  rearrangement 

_  _  <»2N-ln  juit. 

S(u),T)=    NA4  /T    +       E    Aa     E<    E  E       e 

1=2  n=l     m=l 


-  jtut. 
e  >PrCN    =     1]/T 


(D-7) 

and  conditioning  the  interarrival  time  distribution  over  the 

number  o-P  pulses  N  in  the  interval  T 

_  2  <d   1-1   n 

S(u),T)  =  NA4  /T  +  Aa  E    E    E  C0.  (ju)|N=1)  + 

1=2  n=l  m=l 

0." ( jw|N=l) JPrCN  =  13/T 

(D-8) 
where 

jart. 
0.  (  jo)|  N=l  )  =  ECe     3 

is  the  characteristic  -Function  o-F  the  mfc  h  interpulse  spacing 

given  exactly   1  pulses   in  the   interval.      The   -Following 

assumption  is   now  made;   that  the  rnth  interpulse  spacing  is 
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the  sum  of  m  statistically  independent  interpulse  sep- 
arations, and  is  independent  o-P  1,  the  number  of  pulses  in 
the  interval.  This  assumption  will  give  an  asymptotically 
correct  result  for  the  PSD  as  the  observation  interval  T, 
becomes  much  greater  than  p,  the  mean  value  of  the 
interarrival  times. 

With  the  above  assumption 

m 
0.  (j(d|N=1)  =  0  (ju))  (D-9) 

where  0(  ju)    is  the   characteristic  -Function   of  the   mter- 

pulse  separations.   Eqn.  D-8  then  sums  to 

n+1 
_  _  2  <b   1-1   0  ( jcu) -0  ( joj) 

S(o),T)  =  NA4  /T  +  Aa  E    EC — —  + 

1=2  n=l       l-0(j(u) 


n+1 


0-  (  joj)  -0"  (  j(d) 


1  -0"  (  id)) 


UPrCN  =  1]/T 


(D-10) 


and  summing  again 


_  _  2     CD  10  ( jU)) 

S(o),T)     =    NA4  /T    +    Aa     E     < 

1=2 


1  +  1 


0  ( jU)) 


0  (  joj) 


2  2 

(l-0(j(l)))  (l-0(j(l)))  (l-0(j<i))) 


1  +  1 


10"    (   10)) 


0"    (j(D) 


0"    (JO)) 


2  2 

<  1-0"    <  j<i)>  >  (l-0"(jd)))  (l-0"(jU))) 

Rearranging 

_  _  2         _  eo  0  ( ju)) 

S((D,T)     =    NA«  /T    +    Aa  /T-CN    -    PrCN=ll>    E     -C 

1=2 


>PrCN    =     11/T 


(D-l 1 ) 


0"  (  id)) 


(  1"0(  j(D)  )  (  1-0C  j(D)  ) 
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2  0  ( joi)  0"  ( joj) 

-    Aa/T    <1     -    PrCN=0]-PrCN=l>        -C +       > 

2  2 

<  l-0(  jcd)  )  (1-0-  (  jcj)  ) 


1+1  1+1 

2         oo  0  ( joi)  0"  (jut) 

+    Aa  /T    E    PrCN    =     1]        < + >. 

1=2  2  2 

(1-0(joj))  (1-0"  <jw))  (D-12) 

Since  PrCN   =  1]   is  less   than  1   and  |0(u)|   <  1  [Re-f. 

41 :p.   1153    for   any    distribution,   then    the   absolute 

convergence  o-f   the   third   term   is   guaranteed.     For   an 

equilibrium  renewal  process 

N  =  ECN(t) ]  =  T/p  (D-12) 

where  jj   is  the   mean  time  between  renewals  CRe-F.  32:p.  46]. 

Now  letting   T,  the  observation  interval,  go  to  m-Finity  the 

asymptotic  PSD  is 

_  _  AaA4  0(ju)) 

S((u)     =     lim       S(u),T)     =    -C    +  2    ReC     ]     >. 

T  =  >ob  p  2  1-0(jcjj) 

Aa  (D-13) 

This  is   the   expression   for   the   PSD   o-f   an   independent 

increment  point   process   in   terms   o-P   the   characteristic 

•function  o-f  the  interamval  distribution.   A  similar  result 

result  was   developed  by   Cox  and   Miller  by  considering  the 

renewal  intensity  -function  CRe-f.  55]. 
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APPENDIX  E 
SPECTRUM  OF  A  BRANCHING  RENEWAL  POINT  PROCESS 

The  complex   envelope  representation  of  -filtered  impulse 

noise  arising   -Prom  a   branching  renewal  process  on  a  -Finite 

interval  T  is  given  by 

M(T)   N.       -j8.„ 
n(t)  =   E     E  a. „  e      h(t-t.„-T. )         (E-l) 
m=l   n=l 

where  T.    is  the   beginning  of   the  mth  main  interval,  N.  is 

the  number   of  impulses   in  the   mth  interval,   arriving   at 

times  -Ct,  „  >   after  its   onset  and   M(T)  is  the  unit  counting 

process   defining    the   number    of   main    events   in   the 

observation  interval   T.    The  pulse   amplitudes   <a. „>   and 

phases  <9-n>  air&    statistically  independent  random  variables 

and  the   <8.„>  are   taken  to   be  uniform  on  <0,2n>.   The  net 

filtering  effects   on  the   impulses  are  accounted  for  by  the 

time  invariant  complex  envelope  impulse  response,  h(t). 

The  envelope  squared  of  the  process  is 

2  M(T)  N.   2  2 

E  (t)  -  n(t)n-  (t)  =   E     E  a  .  „  |  h  ( t-t.  „  +T.  )  |     (E-2) 

m=l   n=l 

and  using   the  same  assumptions  as  Appendix  D  the  mean  value 

of  the  Bartlett  estimate  of  the  PSD  is 

2 
S(u,T)  =  EC  |  E2  (u),T)  |   D/T  (E-3) 

where  Ea  ( cu  ,  T  )   is  the   length  T   Fourier  Transform  of  E2  (t)  . 

Speci  f ica 1 ly , 
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M(T)        N.     2  -j(j-b.„        -jojT. 

£a  (o),T)=       E  E    a    „  „     e  e  Ha  (a)) 

m=l        n=l 


(E-4) 


with    Ha(cj)     = 


2    -  jort 
|h(t)|     e  dt ,     the    Fourier    trans-Form    oF    the 


0 


magnitude  squared  impulse  response.   To  Find  the  mean  o-P  the 
estimate  we  must  evaluate 

_  2  M(T)  N.  M(T)  N, 

S(M,T)  =  EC  |H,(u)|     E    E    E    E 

m=l  n=l  1=1  k=l 

2    2   -j(u(t.  „-t,  k)   -  ju)(T.  -T>  ) 
a.  „  a,  k   e  e  ]/T 

(E-5) 

In  order   to  simpli-Py   this  expression,   we  make   use  of 

assumptions  1,   2  and   4  -from   section  3D;  (1)  ECa2.„a2i  k  ]  = 

a4  ,  which   will  be   normalized  to  one  -For  simplicity,  (2)  N. 

and  N|    equal  either   Nl  or   N2  alternating   between  the  two 

values  and   (3)  -P(t.„)   and  -P(ti  „  )  ,   the   interarrival   time 

distributions,  alternate   between  the   negative  and  positive 

phase  distribution   parameters.    Furthermore,  let   M  be   an 

even  number   o-P  main   process  points  occurring  at  intervals, 

T„/2,  oP   a  Fundamental   Frequency  (in  our  case  F0  =  uj0  /2n  = 

1/T0  =   60  Hz,   the  power-line   Frequency)  and  let  M2  =  M/2. 

This  assumption   will  greatly   simplify  the  calculations  and 

For  M  greater  than  10,  not  si gn i Ficant ly  reduce  the  accuracy 

oF  the  result. 
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Now   looking    only   at   the   -Factors   involved   in   the 
summation  in  Eqn.  E-5,  six  separate  terms  can  be  identified: 

Term  1   -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  Nl  pulses 

Nl   Nl      jd)(t„  -t„  ) 
=  M2  E    E   ECe  ] 

n=l  k=l 

(E-6a) 

Term  2   -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  N2  pulses 

N2   N2     jui(t„  -t„  ) 
=  M2  E    E   ECe  1 

n=l  k=l 

(E-6b) 

Term  3   -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  Nl  pulses 

M2   M2   -ju)(m-l)T„     Nl   Nl      jwt„      -jut* 
=  E    E   e  E    E   ECe     ]ECe       3 

m=l  1=1  n=l  k=l 

mOl 

(E-6c) 

Term  4  -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  N2  pulses 

M2       M2       -ju)(m-l)T0  N2       N2  ju)t„  -j<i>th 

=    E  E       e  E  E       ECe  DECe  1 

m=l     1=1  n=l     k=l 

mOl 

(E-6d) 

Term  5   -  The  intra-pulse  group  summations  between  pulses  in 

groups  with   Nl  pulses  per  group  to  pulses  in  groups  with  N2 

pu 1 ses 
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jcdCTo /2-TD  )     M2       M2       -ju(m-l)T,        Nl       N2  jut,  -ju>t„ 

=    e  E         E       e  E         E       ECe  3ECe  ] 

m=l     1=1  n=l     k=l 

CE-6e) 

Term  6   -  The  intra-pulse  group  summations  between  pulses  in 

groups  with   N2  pulses  per  group  to  pulses  in  groups  with  Nl 

pulses 

-jw(T0  /2-T*  )  M2   M2   -jw(m-l)T0  N2   Nl      jurt„      -j<i)t„ 
=  e  £         Z       e  E    E   ECe     HECe      ] 

m=l  1=1  n=l  k=l 

(E-6f ) 

It   is    important   to    note   the   expression   -For   the 

characteristic  -Function   of  the   interarrival  times   differs 

between  the   inter-pulse  and   intra-pulse   cases.     In   the 

inter-pulse  case,   discussed  in  Appendix  D,  the  summation  is 

only   over   a   function   of   the   difference   in   the   pulse 

positions  where   in  the   intra-pulse  case   the  summation   is 

over  the  absolute  position  of  each  pulse  in  its  group. 

Terms  1   and  2   can  be   evaluated  using  the  results  from 

Appendix  D  and  assumption  3  from  section  3D: 

Term  1  =  M2  5(w,Nl)  <E-7a) 

Term  2  =  M2  5(<d,N2)  (E-7b) 

and  the  remaining  terms  can  be  determined  by  straightforward 

evaluation  of  the  finite  sums 

cos(M2  2no)/u)0  >  -  1 

Term  3  =  -C  -  M2  >  fi(u,Nl) 

cos  (2no)/u)o  )  -  1  (E-7c) 


cos(M2  2nu)/u)0  )  -  1 

Term  4  =  -C  '■ -  M2  >  Q(id,N2) 

cos(2na)/cD0  )     -     1  (E-7d) 
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Terms  5   and  6  ar&    complex  conjugates  and  can  be  combined  to 

■Form  one  term 

cos(M2  2n_/u>0  )  -  1 

Term  5a  =  2cos(id(T,  /2-T0  )  )  <  — ■ — — — —      >i4/(0,  ,  Nl  ,0a  ,N2) 

cos  (2n_/_0  )  -  1 

(E-7e) 

where  0,   and  0a  ar^&    the  characteristic  Functions  associated 

with  the  interarrival  times  For  each  phase  and  Nl  and  N2  are 

the  number   oF  renewals   associated  with   each  phase.      The 

Functions  5,  Q,  and  iy  are  defined  below; 


2  N 

N  -  N|0(ju)  I  0(ju>)   (1  -  0(j(d)  ) 

5(0, N)  =  _____  _   2  Re-C — — —  >, 

2  2 

I  1  -  0(  j_)  I  (  1  -  0(  joj)  )  (E-8) 


ft(0,N)  = 


N+l  2 
I  0  (  joj)  -  0  (  j_)     J 

2 
I  1  -  0(jo»)  I 


(E-9) 


and 


\+/(0t  ,N1  ,0a  ,N2)  = 


Thu« 


Nl  +  1  N2+1 

(0i  (j_)   -  0,  (j_)      )   (0a  ("J-)   -  0a  ("J-))      ) 


(1-0,  (  JOJ)  )   (1   -  0a  ("J-)  ) 


(E-10) 


S((D,T)  =  -CECHa  (oj)  ]/T>   E  Term(n) 

n=l 


(E-l 1 ) 
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APPENDIX  F 
NON-HOMOGENEOUS  POISSON  PROCESS 

The  -Pol  lowing  derivation  is  for  the  mean  and  variance  of 

a  compound   periodic  non-homogeneous   -Filtered  Poisson  point 

process.     This   proof   is   a   generalization   o-f   a   proof 

presented  in  Papoulis  CRef.  4 1  :  pp .  382-383].   The  expression 

for  such  a  process  on  an  interval  0  to  T  is  given  by 

N(T) 
n(t)  =    Z    at   h(t-tt  ) ,  CF-1 ) 

i  =  l 

where  N(T)   is  a   unit  counting   process  with   periodic  rate 

parameter  X(t) ,   at   is  the  weighting  of  the  ith  point  and  is 

independent  and   identically  distributed  (I  ID)  from  point  to 

point  and  h(t)  is  the  time  invariant  impulse  response  of  the 

f  l  Iter. 

The  time  axis  is  divided  into  consecutive  intervals  Ij 

of    length       acx    where       aoc    =    Oj      -    a4  .  i  .        The    number    of    jumps     in 

the    counting    process     in    the    j*  h      interval     Ij      is    given    by    Arrij   . 

If  aoc   is  sufficiently   small,  then   the  contribution  to  the 

total  noise  process  from  the  jfc h  interval  is 

Arrij 

Anj  *   E  a, h(t-Bj),  CF-2) 

i  =  l 

where  iinj    is   a   Poisson   random   variable   with   its   rate 

parameter  approximately  equal  to 

XCttj  )ao .  (F-3) 
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The    moment    -Function    o-P    the    je  h      interval     is    then 

s    E    a,  h (t-ccj  ) 
i  =  l 
a0j   (s)     *    ECe  1 .  (F-4) 

Conditioning    on    Anij  , 

Amj 
s    E    at  h  (t-<Xj  ) 
i  =  l 
40j  (s)     *     ECe  |Amj     =    k  ]    PrCAm.,     =    k  3  (F-5) 

and 

k   sat  h(t-aj  ) 
a0j  (s)  *  EC  n  e  :  PrCAitij  =  k  H  .  (F-6) 

i  =  l 

Since  the  a,  's  ar^    independent  -Prom  pulse  to  pulse 

k     sah (t-Oj  ) 
ASj  (s)  *   n  ECe  ]  PrCAitij  =  k],  (F-7) 

i  =  l 

where  the  expectation  is  with  respect  to  the  random 
amplitude  a.  Using  the  Poisson  probability  law  and  recog- 
nizing the  series  expansion  For  the  exponential  gives 

sah  (t-cc j  ) 
X(aj  )Acc-CECe  1     -  1> 

A0j  (s)  *   e  .        (F-8) 

Si  nee 

CD 

n(t)  =   E   An, (t) ,  (F-9) 

l  =  -a> 

a  sum  o-P  independent  random  variables,  then 

00 

0„ (s)  =  n  a0, (s)  (F-10) 

1  =-co 
using  the   convolution  property   oF   the   moment   generating 
Function.   The  cumulant  generating  Function  is  defined  as 
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\pn   (s)  =  lnC0„  (s)  ]  , 
where  In  is  the  natural  logarithm  and 


CF-1 1 ) 


H/„  (s)  =  E  ln[A0,  (s)]. 
i  =  -a> 


(F-12) 


As  acc  =>  0,  then 


W«  (s)  = 


p(a) 


sah(t-a ) 
X(a)  C  e         -  13doda.   (F-13) 


Using  a   similar   argument   the   joint   cumulant   generating 
•Function  For  two  random  variables,  n(t,  )  and  n(ta)  is 


\yn  (s,  ,sa  )  = 


p(a) 


Si  ah(t(  -a)  +  saah(ta-oc) 
XCa)  C  e  -  lldoda 


CF-14) 
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APPENDIX  G 
AUTOCORRELATION  OF  CORONA  NOISE 

The  mean  value  of  the  cyclostationary  process  defined  by 

assumptions  1   and  2   in  section   4D  is   given  by  Eqn.  4-12. 

Substituting  Eqn.  4-20  gives 

oo  2 

NT  1/2  -as 

ECe(t) 1    =  ECa3 


NT  1/2  -as 

E   a„  cos  (a),  n  (t-s)  )  <a/n  )    e     ds.  (G-l) 
n=0 


0 


Expanding  the  cosine  and  integrating 

2  2 
-(Do  n 


NT 
ECeCt):  =  (fln74)  E  C  e 

n=0 


4a 


u)o  n 


E 


1 


a„  cos  ( (Do  nt )  + 


2  2 
-(Do  n    k- 1 
(  )    a„  s  l  n  ( (Do  nt )  1 


1/2  k=l   (2k-l ) ! !       2a 
2(an) 

(G-2) 
1/2 
where  CBn/2)      =  ECa]  and  C2k-1)M  =  1  3  5...2k-l  CRef.  54: 

Eqn.  3.897].    The   sine  term   is   in   quadrature   with   the 

dominant  cosine   term  of   Eqn.  G-2   and  represents   a   phase 

shift  in   the  mean   o-f  the   process   relative   to   the   rate 

function.    This  is   due  to   the  filtering   effect  where  the 

impulse  response   of  the  filter  persists,  thus  causing  a  lag 

in  the   mean  value   of  the  process.   For  the  values  of  a  and 

(D0  considered   in  this   work,  the  sine  term  in  Eqn.  G-2  will 
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be  negligible  relative  to  the  cosine  term.   Physically,  this 

means  that   -For  -Pi  Iter   bandwidths   much   greater   than   the 

-fundamental  -Frequency   of  the   corona  noise,  the  phase  shift 

will  be  insignificant. 

Similarly  the   variance  of   the  envelope  when  considered 

as  a  eye  1 ostat i onary  process  is 

2  2 
-o)o  n 


NT  1  1/2  8a 

VarCeCt)]     *     2fl  EC     —     (a/2n)  e  a„  cos  ( tu0  nt )  1  . 

n=0  2 

(G-3) 

Using  Eqn.  4-14  the  covariance  of  the  assumed  process  is 


Cov(t, ,ta )     =    2fl 


Nn 


2  2 

-OC(ti    -S)  -CC(ta   -S) 


Z    an  cos  (ojo  ns  )  (a/n  )  e  e 

!      n=0 
max (ti , t2 ) 

where  after  rearranging  and  completing  the  square 


ds 
(G-4) 


-a  (t,  -ta  )     /2 

Cov(t, ,t8)     =     (2fl(a/n)e  )     x 


NT  -2a (s- (t, +t» ) /2) 

H       an  cos  (  ojo  ns  )  e  ds. 

J     n=0 
max(t, ,ta )  (G-5) 


Lett l ng 


the  covariance  is 


u  =  (s  -  (t,  +  ta)/2), 


(G-6) 


155 


-a  (t,  -ta  )  /2 
CovCt, ,t2 )  =   ((a/R)2fle  )  x 


CO  2 

'   NT  -2au 

JZ  ancos(a)0n(u+(ti  +t2  )/2)e     du 
J  n  =  0 
|-t,  -ta  |  (G-7) 

which  cannot  be  evaluated  in  closed  form. 

The  stationary  autocorrelation  -Function  is  -Pound  by 
using  Eqns.  G-2  and  G-7  and  the  definitions  given  by  Eqns . 
4-15  and  4-19.  The  term  due  to  the  covariance  in  the  auto- 
correlation -Function  is 


2 
-cct  /2 
2fl(a/n)e       /T0 


Tc 


NT 


E  C  a„  cos  (  ojo  n  (2t+t )  ) 
n=0 
0 


-2au 
cos  (  a)o  nu)  e     du 


+  an  si  n  (a)o  n  (t+T /2)  ) 


00  2 

-2au 
sin(o)o  nu)e      duDdt. 


| t |  (G-8) 

The   sinusoidal   terms   will   integrate   to  zero   with   the 

exception  oF   the  a0   cosine  term   and  the  expression  sirnp- 
1 l Fies  to 


2 
-ocr  /2 
2fl(a/n)e        a0 


00  2 

-2au 


du  . 


It  I 


(G-9) 
This  term  will  be  dominated  by  the  First  exponential  and  is 
approximated  by 
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2 
-oct  /2  1/2 

2fle       a0    (a/8n)      .  (G-10) 

The   stationary    autocorrelation   -function    term   due    to 

eye lostat ionary  mean  value  -Function  is 

To 

ECe(t) ]ECe(t+T) ]  dt ,  (G-ll) 


1 

I  0   *■ 


0 


which  using  Eqn.  G-2  is  equal  to 

2  2 
-(Jo  n 


NT  4a  2  2 

2fl     (1/8)       E       e  an     cosuio  nT       +       2fl(l/4)a0     . 

n=l 

(G-12) 

An  approximate  expression  -For  the  stationary  autocorrelation 

■Function    is 


2 
-oct     /2  1/2  2 

R(t)     *     2GCe  a0<oc/8n)  +     (l/4)a0 


2    2 
-o)o     n 

NT  4a  2 

+    (1/8)       E       e  a„     cos  ( u>o  nr  )  3  . 

n=l 

(G-13) 


157 


APPENDIX  H 
LOCALLY  OPTIMUM  ENERGY  RECEIVER 


The  observations   consist  o-f   a  sequence   o-P   N   complex 
samples.   The  hypotheses  are: 
H0  :  X!  =  nt 


versus 


1/2 
Ht  :  x4  =  nt  +  A    st 


i  =  1  ,  .  . . ,N  (H-l ) 

where  nt  are  an  independent,  identically  distributed  (I  ID) 
complex  noise  samples  with  a  joint  inphase  and  quadrature 
density  p ( nc  , na ) .  The  sequence  st  is  a  complex,  zero  mean 
signal  sequence  with  a  known  variance,  <jt  ,  and  A  is  a  real, 
positive  number.  The  1/2  power  o-f  A  is  chosen  to  -facilitate 
the  derivation. 

The  generalized   likelihood  ratio   CRe-P.  6:p.  585,21]  to 
test  the  hypotheses  described  by  Eqn.  H-l  is 


A  ( x  ,  A ) 


1 

N 

=  n 

i  =  l 

V 

L 

1/2  1/2 

p(xc  i  -A   sc  i  ,x8  i  -A   sfl  i  ) 

p(xCi  ,X8  i  ) 


p(Sc  i  ,  S8  i  )  dSc  i  d< 


(H-2) 


Applying  Eqn.  5-8,  and  differentiating  with  respect  to  A, 


N 
tL0  (x,A)     =    E 
1  =  1 


1/2  1/2 

-sc  i  Tip  (Xci-A       Sci   ,Xn  -A       s9i)/2 


1/2 


-OXc  i  p(xc  i    ,X»  t    )  A 
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1/2  1/2 

i9  i  "Dp  (xCi  -A       Sci   |Xg(  -A       Sa  i  )  /2 


]      p(Sc  i    ,  Sa  i    )dSc  .  dss 


1/2 


■ox9  ip(xCi   ,  x9  1  )  A 


A=>0. 
(H-3) 
Letting    A       go    to    0    and    simultaneously    taking    the    expectation 
o-P    sc  1        and    s9  1        requires    L  '  hospital  '  s       rule    to    evaluate    the 
expression     CRe-F.       563.  A-Pter    taking       derivatives    the       test 

statistic    is 


N 
tL0  (x)     =    E 
1  =  1 


2    2  1/2  1/2 

)c  1      "O    p  (xc  i  -A       Sci  ,Xsi  -A       s8  t   )  /2 


■OXc  ,       p(xc  1    ,X8  i    ) 


2  1/2  1/2 

sc  1  s9  1  "Q    p  C  xc  1  -  A       sc  1   ,  xa  1  -  A       s9  1  ) 

TDxc  i  -ox9  1  p(xc  1   ,  x9  i  ) 


2    2  1/2  1/2 

s9  1     "D    p  (xc  1  -A       sc  i   ,  x9  1  -A       s9  i  )  /2 

]    p  (  sc  1  ,  s9  1  )dsc  t  ds9  1 


■ox9  t     p(xC(   ,x,  1  ) 


A  =  >0. 

(H-4) 
In  general,  the  inphase  and  quadrature  noise  components  will 
be  circularly  symmetric  even  though  not  necessarily 
statistically  independent  CRe-P.  41  :p.  133]  such  that 


with 


2  2     1/2 

■F(r)     =    p(  (xc     +    x8      )  ) 


2  2    1/2 

r    =     ( xc     ■+•    x9      ) 


(H-5) 


(H-6) 
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Now  A  may  be  set  to  0  and 


N 
tL0  (r)     =    E 
1  =  1 


2  2  2  2 

sc  i      [r,f  "(r,  )xcl      -    r,  f  (r,  )xc>      +■    r,      H/2 

3 
r,        P(rt  ) 


sc  i  s8  t  C  r,  P  '  '  ( rt  )  xc  ,  x9  i      -    P  '  ( r  t  )  xc  i  x9  t  ] 


rt      f  (rt  ) 


2  2  2  2 

s9  i      [r,f"(r,  )x8,      -    r1-P'(r1)x9i      +    rt      ]/2 

— — — — —     ]p(Sc  i  ,  ss  t  )dsc  i  ds9  i  . 

3 
r,     f(r,  )  (H-7) 

I  f    Set         and    s«  i        are    assumed       "to    have       equal     variances    and    a 

zero    covariance    the    test    statistic    is    now    written    as 


N  2 

tL0  (r)     =         E    (Xi     gt  (rt  ) 
i  =  l 


(H-8) 


where 


P'  '  (r)  P'  (r) 


g(r)     = 


(H-9) 


IP  the  unknown  signal  components  at^e  assumed  to  have  equal 
variances  Prom  sample  to  sample  then  an  equivalent  test 
statistic  is 


N 
t  'L0  (r)  =    E  g(r»  )  . 
i  =  l 


(H-10) 
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